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Abstract 


Laboratory  studies  must  be  scaled,  if  they  are  to  be 
useful  in  predicting  fluid  flow  behaviour  in  petroleum 
reservoirs.  To  this  end,  an  extended  set  of  scaling  groups 
has  been  derived  which  includes  all  of  the  variables  in  the 
immiscible  fluid  displacement  problem.  Of  particular  note  is 
the  result  that  the  geometric  factor  is  unimportant  when 
comparing  a  model  and  prototype,  provided  that  the 
displacement  is  stable. 

To  verify  this  result,  and  to  investigate  the  effect  of 
stability  on  the  immiscible  displacement  problem,  a  series 
of  experiments  was  undertaken.  Distilled  water  was  used  to 
displace  a  Dow  Corning  200  oil  from  tubular,  unconsolidated 
sandpacks.  The  experimental  results  show  that  the  geometric 
factor  is  not  important,  provided  that  the  displacement  is 
either  stable  or  pseudo- stable  (fully  developed  fingering). 
Furthermore,  it  is  demonstrated  that  neither  the  linear  nor 
the  radial  scaling  groups  are  good  correlating  parameters  in 
the  transition  zone  between  stable  and  pseudo-stable 
displacements  ( 1 3 . 56<Ns<900 ) .  Finally,  the  experimental 
results  corroborate  earlier  experimental  work,  delimiting 
stable  and  unstable  displacements,  undertaken  in  this 
laboratory. 
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1.  Introduct ion 


With  the  ever-increasing  demand  for  energy,  it  is  of 
prime  importance  to  employ  and  expand  the  technology 
currently  available  to  the  oil  industry. 

Waterflooding,  the  most  important  secondary  recovery 
process  which  has  been  utilized  for  decades,  is,  as  yet,  not 
fully  understood.  This  incomplete  Knowledge  of  the  process 
may  lead  to  disastrous  results,  if  the  economics  of  the 
project  are  marginal.  Consequently,  the  effect  of  the 
interplay  of  viscous,  capillary  and  gravity  forces  on  the 
recovery  of  oil  from  a  reservoir  must  be  Known. 


1.1  Purpose  of  This  Study 

The  problem  of  conducting  meaningful  laboratory  tests 
on  a  model,  with  the  purpose  of  relating  the  results  to  a 
prototype,  has  been  the  objective  of  many  researchers.  The 
number  of  variables  which  control  the  displacement  process 
is  numerous  and  the  tasK  is  by  no  means  trivial. 

The  purpose  of  this  study  was  twofold.  The  first 
objective  was  to  establish  a  set  of  scaling  groups  which 
included  all  of  the  variables  in  the  problem  as  it  is 
understood  today.  The  approach  used  to  reach  this  first 
objective  was: 

1)  to  review  the  results  of  the  methods  of 

dimensional  and  inspectional  analysis  found  in 
the  literature, 
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2)  to  derive  the  scaling  groups  in  a 

mul ti -dimensional  system  using  a  current 

capillary  pressure  model,  and 

3)  to  verify  certain  scaling  groups  experimentally. 

The  second  objective  was  to  investigate  further,  and 

verify,  the  criterion  which  determines  the  onset  of 
floodfront  perturbations.  The  specific  areas  of  interest 
were : 

1)  Water-wet,  residual  water  saturation  system. 

a)  Better  locate  the  onset  of  instability 
experimental ly , 

b)  Investigate  the  possibility  of  a  second 
region  in  the  recovery  curve,  at  a  high 
value  of  the  stability  number,  where  the 
recovery  is  insensitive  to  frontal 
perturbations,  and 

c)  Extend  the  data  collected  by  Peters  (1979) 
to  small  and  large  values  of  the  stability 
number . 

2)  Oil-wet,  no  residual  water  saturation  system. 

a)  Attempt  to  locate  the  onset  of  instability 
in  the  recovery  curve,  as  suggested  by  the 
theory. 

b)  Better  define  the  second  region  in  the 
recovery  curve  for  high  values  of  the 
stability  number  and  determine  its 
magni tude . 
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c)  Extend  the  data  collected  by  Peters  (1979) 
to  both  small  and  large  values  of  the 
stabi 1 i ty  number . 

The  purpose  of  extending  the  data  collected  by 
Peters  (1979),  to  small  and  large  values  of  the  stability 
number,  was  to  verify  the  suggestion  that  the  recovery  is 
independent  of  the  stability  number  for  these  magnitudes.  It 
was  also  intended  to  investigate  further  any  possible 
anomalies  which  might  appear  in  the  data. 


1.2  Stabilization  and  Stability-Terminology  Defined 

Throughout  the  literature,  the  terms  'stabilization' 
and  'stability'  appear  often.  In  view  of  the  fact  that  both 
terms  apply  to  the  behaviour  of  the  floodfront,  it  is 
important  to  understand  the  physical  interpretation  of  these 
terms  in  order  to  avoid  confusion. 

'Stabilization'  occurs  when  the  frontal  region  of  the 
saturation  profile  has  attained  a  fixed  configuration 
invariant  with  time.  This  saturation  profile  should  approach 
that  predicted  by  Buckley-Leveret t  (1941)  theory. 

During  the  life  of  a  waterflood,  frontal  per turbat ions 
(viscous  fingering)  may  arise.  The  displacement  is 
considered  to  be  'stable'  if  these  frontal  perturbations  are 
attenuated.  Thus  a  'stable'  displacement  may  be  either 
'stabilized'  or  'unstabilized',  depending  on  the  magnitude 
of  certain  variables  which  control  the  displacement  process. 
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2.  Theory  and  Literature  Review 


2.1  Theory  of  Immiscible  Displacement 

2.1.1  Fundamental  Equations 

The  partial  differential  equation  describing  the 
displacement  of  one  fluid  by  another  through  a  porous  medium 
may  be  derived  by  combining  Darcy's  law  and  the  equations  of 
continuity  for  each  phase. 

Darcy's  law,  in  a  general  three-dimensional  form,  for 
the  simultaneous  flow  of  two  immiscible  and  incompressible 
fluids  through  a  homogeneous  and  isotropic  porous  medium  may 
be  wri tten  as : 


W 


(arad  P  , 
w 


2.  1 


and 


MV 


y 


o 


(grad  PQ 


2.2 


for  each  phase. 

If  the  porous  medium  is  considered  to  be 
non-deformable,  then  the  equations  of  continuity  may  be 
written  as  (Buckley  and  Leverett,  1942): 


di v  v  = 
w 


as 


-  <j> 


w 


at 


2.3 


and 


4 
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as  as 

o 


2.4 


for  each  phase.  Note  that  conservation  of  mass  requires  that 
the  saturations  must  sum  to  unity,  viz: 


2.5 


A  capillary  pressure  definition  (Leverett,  1941) 


2.6 


P  (S  )  =  P  -  P 
C  W  0  w 


is  required  to  link  Equations  2.1  and  2.2. 

2.1.2  Relative  Permeability  and  Mobility  Ratio 

Before  proceeding  with  the  derivation  of  the  partial 
differential  equation,  it  is  convenient  to  introduce  the 
concepts  of  relative  permeability  and  mobility  ratio. 
Relative  permeability  is  defined  as  the  effective 
permeability  to  a  specific  fluid  divided  by  some  base 
permeability,  viz: 


2.7 


and 


2.8 
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There  are  three  different  base  permeabilities  which  are 
generally  used  throughout  the  literature  (Craig,  1971): 

1)  The  absolute  air  permeability; 

2)  The  absolute  water  permeability; 

3)  The  permeability  to  oil  at  residual  water 
saturation . 

Equations  2.1  and  2.2  both  contain  a  term  of  the 
form  (k/y  ) .  This  term  is  defined  as  the  mobility  of  the 
fluid.  As  will  be  demonstrated  in  a  later  section,  a 
mobility  ratio  will  be  defined.  Suffice  it  to  state  the 
general  form  of  the  ratio  of  mobilities  as: 


W 


u 


w 


uo 

W 


2.9 


The  ratio  of  mobilities  may  be  defined  as  the  water  mobility 
divided  by  the  oil  mobility.  Using  the  definitions  of 
relative  permeability  given  by  Equations  2.7  and  2.8  yields: 


M 


P. 


Kwb 

TOT  Kob 


2.10 


A  general  mobility  ratio  may  be  defined  as: 


M 


g 


^wb  yo 
^ob  uw 


2.11 


It  has  been  stated  in  the  literature  (Bentsen,  1976) 
that  defining  the  relative  permeability  in  terms  of  a  single 
base  permeability  may  not  be  desirable.  If  this  is  done, 
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then  a  non-dimensional izing  procedure  will  select  the 
viscosity  ratio  as  the  scaling  group.  However,  it  is 
generally  thought  that  the  mobility  ratio  is  a  better 
scaling  group  (Craig,  1971). 

After  the  manner  of  Perkins  and  Collins  (1960),  in 
which  a  different  base  permeability  is  used  to  normalize  the 
two  effective  permeabilities,  Equations  2.7  and  2.8  may  be 
written  as: 


2.12 


and 


2.13 


where  K.„  and  K  are  the  permeabilities  to  water  at 
residual  oil  saturation  and  to  oil  at  residual  water 
saturation  respectively.  Furthermore,  these  relative 
permeabilities  may  be  defined  in  terms  of  a  normalized  water 
saturation  (  S*  )  : 


2.14 


K  (S*) 
K 


2.15 


Figures  2.1  and  2.2  illustrate  typical  and  normalized 
relative  permeabilities,  respectively. 
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FIG:  2. 1  TYPICAL  RELATIVE 

PERMEABILITY  CURVES 


FIG.  2.2:  TYPICAL  NORMALIZED 
RELATIVE  PERMEABILITY  CURVES 
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2.1.3  Partial  Differential  Equation 

The  equations  describing  the  fluid  displacement  problem 
may  now  be  combined  (see  Appendix  A.1  for  details)  to  yield 
the  result: 
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2.  17 


and 


Ap  =  pw  -  p0  2.18 

This  result  is  a  generalization  of  the  non-linear, 
parabolic  partial  differential  equation  derived  by 
Rapoport  (1955)  and  is  not  subject  to  an  analytical 
solution.  The  equation  accounts  explicitly  for  the  viscous, 
gravitational,  and  capillary  forces. 

2.1.4  Fractional  Flow  Equation 


If  the  fractional  flow  of  water  is  defined  as: 
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2.19 


then  Equation  A. 1.15  becomes: 
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2.20 


and  Equation  2.3  becomes: 


div  f 
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2.21 


Equations  2.20  and  2.21  also  describe  the  fluid 
displacement  problem. 

2.1.5  Capillary  Pressure  Models 

Several  capillary  pressure  models  have  been  proposed 
throughout  the  literature.  Of  particular  note  are  the: 


1)  Leverett  J-function  (1941) 


2.22 


2)  Modified  Leverett  J-function  (Rose  and 
Bruce,  1949) 


2.23 
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3)  Bentsen  and  Anli  (1976) 
PC(S*)  =  -  alnS*  +  Pd 


2.24 


The  J-function  attempts  to  correlate  the  physical 
properties  of  the  rock  and  fluid.  The  expression  given  by 
Equation  2.23  includes  the  wetting  angle.  The  J-function  is 
restricted  to  correlating  satisfactorily  data  obtained  from 
unconsolidated  porous  media  or  certain  lithological 
classifications  (Brown,  1951).  Bentsen  and  Anli  (1976) 
proposed  the  model  given  by  Equation  2.24  in  an  attempt  to 
overcome  the  restrictions  of  the  J-function. 

The  capillary  pressure  normalizing  parameter 
incorporates  the  effects  of  interfacial  tension,  wettability 
and  pore-size  distribution.  Both  the  displacement  pressure 
and  capillary  pressure  normalizing  parameter  must  be 
determined  experimentally. 

In  normalized  form,  the  capillary  pressure  equation  may 
be  wri tten  as : 


IT 
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(S*)  =  -  1 nS* 
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where 


2.26 


or 
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pc($*)  =  <™C(S*)  +  Pd  2.27 

Figures  2.3  and  2.4  illustrate  typical  and  normalized 
capillary  pressure  curves,  respectively.  In  physical  terms, 
the  area  under  the  capillary  pressure  curve  represents  the 
work  done  in  creating  new  boundary  surface. 

In  view  of  the  fact  that  the  model  parameters  are 
available  experimentally,  the  Bentsen  and  Anli  (1976)  model 
will  be  used  throughout  this  study  in  any  mathematical 
derivations  presented.  Moreover,  Golaz  (1979)  showed  that 
the  capillary  pressure  normalizing  parameter  retained  the 
same  value,  regardless  of  the  capillary  pressure  model  used. 


2.2  Immiscible  Displacement  Mechanisms 

The  fluid  displacement  problem  introduced  in  Section 

2.1  may  be  represented  by  various  displacement  mechanisms. 
The  Muskat  model  is  the  simplest  and  easiest  to  apply  in  an 
analysis.  The  mathematical  development  of  the  stability 
theories  is  based  on  the  Muskat  theory. 

2.2.1  Muskat  Displacement  Mechanism 

With  the  assumption  that  the  region  in  which  an 
appreciable  saturation  gradient  exists  is  small,  then  the 
porous  medium  can  be  divided  into  two  parts:  one  containing 
only  oil  and  irreducible  water,  and  the  second  containing 
only  water  and  immobile  residual  oil.  Furthermore  in  the 
first  region,  which  is  located  ahead  of  the  front,  only  oil 


. 

' 


13 


FIG:  2.3  TYPICRL  CRPILLRRY 

PRESSURE  CURVE 


FIG.  2.4:  TYPICRL  NORMALIZED 
CRPILLRRY  PRESSURE  CURVE 
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is  mobile;  and  in  the  second  region,  which  is  located  behind 
the  front,  only  water  is  mobile.  Figure  2.5  illustrates 
schematically  the  saturation  distribution  in  such  a  system. 

The  surface  separating  these  two  regions  may  be 
interpreted  as  a  saturation  discontinuity.  The  velocity  of 
this  interface  is  given  by: 


v 
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vwn 

♦t1  -  Sor  *  Vi> 


2.28 


Applying  Darcy's  law  and  the  equation  of  continuity,  and 
providing  that  gravity  and  capillary  pressure  is  neglected, 
the  problem  may  be  formulated  as: 
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FIG.  2.5:  SATURATION  DISTRIBUTION 
FOR  THE  MUSKAT  AND  BUCKLEY-LEVERETT  MODELS 
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Integrating  Equations  2.29  and  2.30,  and  applying  the 
boundary  conditions  result  in: 

pw  =  Cxx  +  C2  2.35 

p0  =  CJx  +  C£  2.36 


where 


Cl  '  -  M'r'L  +  (f  -  hr)xf  2'37 

C2  =  Pj  2.38 

M  aP 

c>  =  '  MrL  +  (1  -  Mr}xf 
(1  -  M  )xfaP 

C2  =  MrL  +  (1  -  Mr)xf  +  Pl 
aP  =  P!  -  P2  2.41 


2.39 

2.40 


The  rate  of  advance  of  the  front  may  be  readily  derived 
from  these  equations  and 
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Hence 
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dt 
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Defining 
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Equation  2.44  may  be  integrated  wi th  Kf  =  0  @  t*  =  0  to 
yield : 
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-  VMr  -  2MMr  “  1)t* 
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2.47 


Figure  2.6  illustrates  the  fractional  distance 
traversed  by  the  front  for  a  linear  displacement  showing  the 
effect  of  mobility  ratio  (Muskat,  1937). 


2.2.2  Buck  ley- Leveret t  and  Dietz  Displacement  Mechanisms 

Under  certain  circumstances,  the  piston- like 
displacement  discussed  in  the  previous  section  may  not 
describe  the  process  adequately.  The  determination  of  the 
water  saturation  at  breakthrough  is  fundamental  to  the 
establishment  of  over-all  reservoir  production.  If  the 
saturation  in  the  trailing  zone  is  not  immediately  reduced 
to  the  residual  oil  saturation,  then  the  Muskat  model  cannot 
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be  applied.  This  trailing  zone  is  illustrated  in  Figure  2.5. 

There  are  two  existing  theories  of  immiscible 
displacement  that  have  found  acceptance  in  the  literature. 
According  to  the  theory  suggested  by  Buckley  and  Leveret t 
(1942),  the  invading  water  progresses  with  the  formation  of 
a  front  and  a  trailing  zone  (or  transition  zone  with  an 
appreciable  saturation  gradient).  With  the  assumption  that 
gravity,  capillarity  and  variations  in  density  are 
negligible,  Buckley  and  Leverett  (1942)  wrote  the  following 
equation  for  the  one-dimensional  flow  of  two  immiscible 
fluids: 


<j)A  9X 


2.48 


The  above  equation  may  be  derived  from  a  material  balance 
and  Darcy's  law.  From  Equation  2.48,  the  rate  of  advance  of 
the  front  may  be  formulated  as: 


dx 

dt 


df 
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2.49 


Equation  2.49  is  the  so  called  '  Buckley-Leveret t'  equation. 
The  theory  assumes  that  the  oil  and  water  flow 
simultaneously  through  the  porous  medium.  Furthermore,  the 
relative  permeability  curves  must  be  known  explicitly. 

The  theory  of  Dietz  (1953)  is  designed  for  a 
two-dimensional  flow  system  and  assumes  that  the  water 
invades  the  formation  in  the  form  of  a  tongue.  It  is  further 
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assumed  that  the  oil  saturation  in  the  invaded  zone  is 
instantaneously  reduced  to  the  residual  oil  saturation. 
Unlike  the  BucKl ey- Leveret t  theory,  the  oil  and  water  flow 
separately:  the  oil,  ahead  of  the  front  through  a  porous 
medium  of  constant  relative  permeability  to  oil;  and  the 
water,  in  the  invaded  zone  with  a  constant  relative 
permeability  to  water  at  residual  oil  saturation.  The 
saturation  distribution  for  the  Buckley-Leveret t  and  Dietz 
models  is  illustrated  in  Figure  2.7. 

The  analogous  expressions  to  Equations  2.48  and  2.49 
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Y  =  y/h 
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=  coordinate  of  the  oil/water  interface  as  a 
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fraction  of  formation  thickness;  its  physical 
meaning  being  restricted  to  0<y<  h 
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2.54 


=  the  difference  in  saturation  between  the  oil  zone 
and  water  zone 


The  Dietz  theory  is  simpler  to  apply  than  the 
Buckley- Leveret t  theory  in  view  of  that  fact  that  only  the 
end-point  relative  permeabilities  need  be  known.  Because  the 
Dietz  theory  is  not  restricted  to  assuming  a  single  tongue, 
the  differences  in  the  two  theories  diminish  if  a  great  many 
tongues  exist.  The  two  theories  are  identical,  provided  the 
relative  permeabilities  are  defined  as  linear  functions  of 
saturation  (  =  S* ,  Kro  =  1-S* )  . 


2.3  Scaled  Models  and  Scaling  Groups 

Dimensional  models  were  first  introduced  to  reservoir 
engineering  by  Leverett  and  co-workers  (1942)  in  order  to 
study  the  water-drive  process.  The  resulting  set  of  scaling 
groups  enables  the  experimenter  to: 

1)  vary  the  pertinent  variables  of  the  problem  in  a 
systematic  way  as  a  means  of  evaluating  the 
process  and 

2)  relate  laboratory  results  to  field  applications. 
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2.3.1  Derivation  of  Scaling  Groups 

The  scaling  groups  may  be  derived  by  either  dimensional 
analysis  or  inspect ional  analysis  techniques. 

The  dimensional  analysis  approach  consists  of 
considering  all  the  variables  involved  in  the  process  being 
studied,  and  combining  these  variables  into  dimensionless 
groups  as  dictated  by  Buckingham's  (1915)  pi -theorem.  In  an 
inspectional  analysis,  all  equations  which  describe  the 
process  of  interest  are  combined  to  form  a  single  equation. 
The  parameters  form  the  scaling  groups.  If  any  of  the 
equations  are  of  the  differential  type,  then  inclusion  of 
initial  and  boundary  conditions  is  also  required. 

The  dimensionless  groups  resulting  from  such  an 
analysis  may  be  divided  into  the  following  classes: 

1)  independent  groups,  or  the  dimensionless  form  of 
the  independent  variables; 

2)  dependent  groups,  or  the  dimensionless  form  of 
the  dependent  variables;  and 

3)  similarity  groups,  which  are  independent 
constant  groups. 

Dimensional  analysis  does  not  require  that  the  process 
be  expressed  by  equations.  Hence,  complete  knowledge  of  the 
process  is  not  required.  This  is  a  distinct  advantage  over 
inspectional  analysis.  Special  care  must  be  taken  to  include 
all  of  the  pertinent  variables  associated  with  the  problem 
at  hand.  Furthermore,  the  resulting  dimensions  less  groups 
may  not  have  any  clear  physical  interpretation  and,  in 


. 
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general,  the  groups  are  not  unique.  An  inspectional  analysis 
of  all  the  equations  usually  results  in  the  scaling  groups 
having  a  clear  physical  meaning. 

The  early  work  of  Leverett,  et  al.  (1942)  has  been 
built  upon  and  used  extensively.  Presented  in  Appendix  A. 2 
are  the  details  of  the  derivation  of  the  scaling  groups  for 
the  immiscible  fluid  displacement  problem,  using  the  two 
methods  mentioned  above. 

Engleber ts  and  Klinkenberg  (1951)  extended  the 
dimensional  analysis  of  Leverett,  et  al.  (1942).  Their 
analysis  yielded  fewer  pi -terms  and,  coupled  with  simpler 
experimental  systems,  obtained  results  of  more  general 
signi f icance. 

Rapoport  and  Leas  (1953)  presented  their  experimental 
results  for  linear  waterfloods  in  which  they  systematically 
varied  length,  wetting  phase  viscosity,  and  flow  rate. 

Rapoport  (1955)  was  the  first  to  present  the  derivation 
of  the  scaling  groups,  using  inspectional  analysis,  for  a 
three-dimensional  cartesian  coodinate  system.  Geertsma,  et 
al .  (1956)  derived  the  dimensionless  groups  for  cold-water 
drive,  hot-water  drive  and  solvent  injection  by  inspectional 
analysis.  They  then  'completed'  the  resulting  set  of 
pi -terms  by  means  of  dimensional  analysis.  This  was  done  in 
an  attempt  to  combine  the  advantages  of  both  techniques. 

Perkins  and  Collins  (1960)  recognized  that  previous 
derivations  of  the  scaling  groups  required  that  the 
functional  relationships  of  both  relative  permeability  and 
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capillary  pressure  be  identical  in  the  model  and  prototype. 
They  presented  a  new  set  of  scaling  groups  in  which  the 
functional  relationships  were  allowed  to  be  different.  This 
was  done  by  redefining  relative  permeability  and  normalizing 
the  water  saturation.  Consequently,  the  new  set  of  scaling 
groups  demanded  that  the  functional  dependence  of  relative 
permeability  and  capillary  pressure  be  on  normalized  water 
saturation  instead  of  water  saturation. 

Bentsen  (1976)  improved  the  work  of  Perkins  and  Collins 
(1960)  for  a  linear  displacement  problem  using  the  method  of 
Heliums  and  Churchill  (1961).  The  method  used  is  essentially 
a  modified  inspect ional  analysis.  However,  the  resulting 
scaling  groups  were  difficult  to  interpret  physically. 
Consequently,  Bentsen  (1976)  redefined  the  pi -terms,  thus 
attaching  a  physical  significance  to  the  dimensionless 
groups.  A  derivation  extending  the  work  of  Bentsen  (1976)  to 
a  cylindrical  coordinate  system  is  offered  in  the  next 
section . 


2.3.2  Dimensionless  Form  of  the  Fluid  Displacement  Equations 

As  a  consequence  of  applying  a  modified  inspect ional 
analysis,  Equations  2.20  and  2.21  take  on  the  dimensionless 
form  ( Append i x  A . 2 . 3 ) : 
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In  light  of  the  previous  discussion,  the  following 
classes  of  dimensionless  groups  arise: 

1 )  Independent  groups 

e,  t 


2)  Dependent  groups 


3)  Similarity  groups 


In  addition  to  these  groups  there  are  the  functions 
Kpv/S*)  »  KrQ(S*)  and  ^(S*)  •  These  have  the  properties  of 
similarity  groups  and,  consequently,  must  be  the  same  in  the 
model  and  prototype. 

If  the  displacement  is  essentially  linear,  then  the 
saturation  gradient  in  the  radial  direction  is  small.  Hence: 
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and 
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Furthermore,  the  pressure  at  any  given  cross-section  will  be 
uniform.  This  implies  that: 
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Hence  Equations  2.55  and  2.56  reduce  to: 
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and 


fw  =  G(s*)  -  Ncc(s*)  If  2-74 

The  above  two  equations  are  identical  to  those  derived  by 
Bentsen  (1976)  and  apply  only  if  there  is  no  transverse 
flow. 

2.3.3  Remarks  on  Scaling  Groups 

It  is  evident  that  in  the  design  of  a  model,  it  is  a 
difficult,  if  not  impossible,  task  to  meet  the  requirements 
of  all  the  similarity  groups  simultaneously.  It  is  desirable 
then  to  relax  the  groups. 

Perkins  and  Collins  (1960)  showed  that  redefining  the 
relative  permeability  curves  in  terms  of  the  normalized 
water  saturation  was  less  restrictive  than  previous 
definitions.  Van  Daalen  and  van  Domeselaar  (1972)  presented 
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experimental  evidence  indicating  that  for  small  transition 
zones  the  functions  K^S*)  ,  KrQ(S*)  and  ^(S*)  did  not 
influence  the  production  behaviour.  Under  these 
circumstances,  they  do  not  have  to  be  satisfied. 

Bentsen  (1976)  suggested  that  the  capillary  number,  if 
sufficiently  small,  may  also  be  neglected.  This  group  may  be 
interpreted  as  a  measure  of  the  relative  importance  of  the 
capillary  forces  to  viscous  forces.  It  was  demonstrated  both 
theoretically  (Bentsen,  1978)  and  experimentally  (Saeedi, 
1979)  that  for  a  capillary  number  of  less  than  0.01,  the 
recovery  was  independent  of  this  group.  The  dependence  of 
the  recovery  on  the  capillary  number  is  illustrated  in 
Figure  2.8,  and  the  effect  of  mobility  ratio  is  also  shown. 

The  ratio  aKwr/|v|ywD  ,  similar  to  Nc£  but 
representative  of  the  capillary  to  viscous,  forces  in  the 
radial  direction.  This  ratio  may  be  neglected  without 
restriction  where  no  crossflow  occurs.  However,  when  gravity 
and  viscous  forces  play  an  important  role  (gravity  underride 
and  viscous  fingering),  crossflow  is  present.  This  group 
must  then  be  taken  into  account.  If  the  magnitude  of  the 
radial  capillary  number  is  small,  say  of  the  order  of  0.01, 
then  it  may  be  anticipated  that  the  fingers  will  have 
stabilized  in  the  same  sense  that  the  floodfront  is  said  to 
be  stabilized  (in  the  direction  of  bulk  flow)  when  the 
linear  capillary  number  is  small. 

Finally,  it  should  be  noted  that  the  recovery  at 
breakthrough  for  the  immiscible  fluid  displacement  problem 
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is  some  complicated  function  of  mobility  ratio,  gravity 
number,  capillary  number,  geometric  factor,  and  the  effects 
of  instability  of  the  floodfront.  The  effects  of  instability 
of  the  floodfront  are  discussed  in  a  later  section.  The 
recovery  function  may  be  represented  as: 

Recovery  =  f  (Mr,  Ng,  Nc£,  L/D,  N$) 
where  N$  accounts  for  the  effects  of  floodfront  instability. 
If  the  displacement  is  stable  and  stabilized,  then: 

Recovery  =  f  (M^,  N  ) 

Furthermore,  if  gravity  can  be  neglected,  then  the  simplest 
case  is: 

Recovery  =  f  (M  ) 

2.4  Stability  Theory  of  Immiscible  Displacement 

The  theory  presented  in  the  preceding  sections 
implicitly  assumes  that  the  floodfront  has  attained  a 
stabilized  configuration  and  that  no  frontal  per turbat ions 
occur.  Consequently,  the  stable  model  is  incapable  of 
predicting  the  waterflood  performance  of  a  prototype  in  a 
regime  where  viscous  fingering  dominates  the  recovery 
process.  As  a  result,  the  instabilities  associated  with  the 
model  may  not  be  properly  accounted  for  in  the  prototype. 
Hence  the  need  arises  to  develop  a  mathematical  model  which 
accounts  explicitly  for  the  floodfront  perturbations. 


. 


32 


2.,-4 . 1 _ Experimental  Observations  and  Mathematical 

Descriptions  of  Instability 

When  one  fluid  is  displaced  by  another  fluid  with  a 
lower  viscosity,  the  fluid-fluid  interface  is  in  an  unstable 
state  and  has  a  tendency  to  'break  up'  .  The  nature  of  this 
instability  phenomenon  was  first  observed  by  Engleberts  and 
KlinKenberg  (1951).  The  streamers  which  appeared  in  the 
visual  models  prompted  them  to  coin  the  term  'viscous 
fingers' . 

For  an  adverse  viscosity  ratio  of  24,  Engleberts  and 
KlinKenberg  (1951)  observed  a  continuing  decrease  in 
recovery  at  breakthrough  for  both  hydrophilic  and  oleophilic 
powders.  This  decrease  in  recovery  was  attributed  to  the 
onset  of  perturbations  at  the  floodfront  and  growth  of  the 
incipient  fingers  toward  the  outlet,  thus  bypassing  oil  in 
the  model . 

Van  Meurs  (1957)  described  a  technique  which  permitted 
complete  visual  observation  of  the  displacement  process.  It 
was  seen  that,  for  a  viscosity  ratio  of  unity,  the 
fluid-fluid  interface  was  stable  and  the  displacement 
efficiency  was  high.  Further  tests  at  an  adverse  viscosity 
ratio  of  80  confirmed  the  observations  of  Engleberts  and 
Klinkenberg  (1951).  Furthermore,  van  Meurs  (1957)  observed 
that  the  production  of  oil  after  breakthrough  is  due  to  the 
growth  of  the  fingers  already  formed.  These  tests  were 
perforned  in  an  oil-wet  system  containing  no  residual  water 
saturation. 
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The  results  of  the  visual  observations  of  the 
displacement  process  prompted  van  Meurs  and  van  der 
Poel  (1958)  to  attempt  to  describe  the  displacement  process 
involving  viscous  fingering  mathematically.  This  model, 
which  parallels  the  theory  of  Buckley  and  Leverett  (1942)  in 
some  respects,  attempts  to  account  for  frontal 
perturbations.  However,  little  recognition  has  been  given  to 
this  theory  in  the  literature. 

Chuoke,  et  al.  (1959)  were  the  first  to  present  a 
theoretical  description  of  viscous  fingering.  With  the 
assumption  that  the  frontal  disturbances  may  be  represented 
by  a  Fourier  decomposition,  and  using  first  order 
perturbation  theory,  they  defined  a  stability  index  as: 


For  an  excellent  detailed  derivation  of  Equation  2.75,  the 
reader  is  referred  to  the  work  of  Peters  (1979). 

In  view  of  the  fact  that  the  sign  of  the  stability 
index  determines  the  stability  classification,  the  necessary 
and  sufficient  condition  for  instability  (n  positive)  is: 

jv  +  (p0  -  Pw)gcosa'  >  Y*S2  2.76 

A  critical  velocity  may  be  defined  such  that  the  left  hand 
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side  of  Equation  2.76  is  equal  to  zero,  that  is: 


/yo  %  \ 

\  or  '  V 


Vc  +  (po  “  Pw)9C0Sa'  =  0 


2.77 


or 


K  Apgcosa' 

\/  =  wr _ _ 

<=  VMr  -  D 

From  Equations  2.76  and  2.77,  there  results: 


(V  -  Vc)  >  Y*S2 


2.78 


2.79 


or 


X  >  2tt 


KwrV* 

-  D(V  -  V  'J 


1/2 


2.80 


which  is  the  condition  for  instability. 

From  Equation  2.80,  a  critical  wavelength  may  be 
defined  as: 


x 


c 


2tt 


-  1KV  -  vc) 


2.81 


Differentiation  of  Equation  2.75  with  respect  to  6, 


and 


equating  to  zero,  determines  the  wavelength  of  the  maximum 
instability,  viz: 
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l\  -W1HV  -  Vc) 


1/2 


2.82 


X  -  2tt/3 
m 


or 


2.83 


To  verify  these  theoretical  findings,  Chuoke,  et  al.  (1959) 
used  three  types  of  experimental  systems:  a  Hele-Shaw  model 
inclined  at  45°,  a  van  Meurs  model  with  no  residual  water 
and  a  van  Meurs  model  with  residual  water. 

Because  in  a  Hele-Shaw  model  y*  is  equal  to  y,  the 
wavelengths  of  the  maximum  instability  can  be  directly 
calculated.  Excellent  agreement  between  the  calculated  and 
measured  wavelengths  was  obtained  in  the  paral lei -plate 
tests . 

The  displacement  tests  showed  two  trends  consistent 
with  the  theory.  It  was  observed  by  Chuoke,  et  al.  (1959) 
that  small  fingers  are  formed,  when  the  oil  viscosity  is 
high.  Furthermore,  the  lower  the  bulk  interfacial  tension, 
the  smaller  the  fingers. 

Chuoke,  et  al .  (1959)  assumed  a  linear  relationship 
between  the  effective  interfacial  tension,  which  is  an 
unknown  in  porous  media,  and  the  bulk  inter  facial  tension  to 


be : 


y*  =  Cy 


2.84 


Hence,  Equation  2.82  becomes: 
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2.85 


where  C  is  Chuoke' s  constant  which  is  of  the  order  of  2tt/3  . 
Chuoke,  et  al.  (1959)  found  that  the  experimental  data 
fitted  Equation  2.85  well  when  C  was  taken  to  be  equal  to 
30.  However,  this  value  of  Chuoke' s  constant  was  only  valid 
for  an  oil -wet  system  with  no  residual  water. 

Insufficient  data  were  accumulated  by  Chuoke,  et 
al .  (1959)  to  arrive  at  a  quantitative  estimate  of  the 
system  constant  in  a  water -wet  porous  medium.  However,  it 
was  speculated  that  Chuoke' s  constant  would  be  much  larger 
in  this  type  of  system. 

In  terms  of  field  applications  when  instability  occurs, 
the  interpretation  of  model  tests  may  be  misleading. 
Consequently,  Chuoke,  et  al .  (1959)  suggested  that  the 

p 

similarity  groups  should  contain  the  additional  group  um/D) 
to  account  for  the  floodfront  perturbations. 

It  was  pointed  out  in  Section  2.3.3  that  it  may  not  be 
feasible  to  satisfy  simultaneously  all  the  similarity 
groups.  Chuoke,  et  al.  (1959)  suggested  that  if  the 
wavelength  is  of  the  order  of  the  diameter  of  the  sandpack, 
exact  scaling  is  required;  otherwise  it  may  be  neglected. 

De  Haan  (1959)  presented  further  evidence  confirming 
the  theory  of  Chuoke,  et  al.  (1959).  In  an  oil -wet  system, 
it  was  observed  that  a  decrease  in  recovery  occurred  with  an 
increase  in  the  scaling  factor,  once  fingering  commenced.  At 
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still  higher  rates,  de  Haan  observed  a  levelling  in  the 
breakthrough  recovery  curve. 

Outmans  (1962)  attempted  an  improvement  to  the 
first-order  theory  by  including  the  non-linear  terms  in  the 
equations  describing  the  conditions  at  the  interface.  Of 
particular  note  is  his  conclusion  that  the  shape  and  growth 
of  the  fingers  cannot  be  calculated  from  the  equations  which 
arise  from  first-order  theory.  Consequently,  the  first  order 
equations  must  be  replaced  by  more  complicated  expressions. 
He  suggests,  if  interfacial  tension  is  negligible,  that  the 
conclusion  from  linear  theory--that  the  rate  of  growth  of 
the  incipient  finger  will  continue  without  bound--is 
incorrect.  Finally,  Outmans  (1962)  questioned  the  validity 
of  the  dynamic  boundary  condition  at  the  interface  suggested 
by  Chuoke,  et  al.  (1959).  This  boundary  condition,  which 
allows  for  capillarity  at  the  interface,  was  redefined  by 
Outmans  ( 1 962 ) . 

Using  a  first-order  analysis,  Rachford  (1964)  presented 
an  extension  to  the  theory  proposed  by  Chuoke,  et  al. 

(1959).  By  including  a  saturation  transition  zone 
( Buckley- Leveret t  type  displacement),  Rachford  (1964) 
concluded  that  laboratory  tests  scaled  by  the  usual  groups 
are  also  correctly  scaled  for  the  presence  of  viscous 
fingering  in  a  water-wet  system. 

This  first  attempt  to  study  the  stability  problem  using 
Buck ley- Leveret t  theory  gave  rise  to  highly  non-linear 
perturbation  equations.  Rachford  (1964),  unable  to  solve 
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these  equations,  violated  the  nonlinearity  of  the 
perturbation  equations  and  consequently  his  results  may  be 
questionable. 

Perkins  and  Johnston  (1969)  studied  immiscible 
fingering  in  linear  Hele-Shaw  and  bead-packed  models.  It  was 
observed  that  viscous  fingering  was  severe  in  the  system 
containing  no  residual  water,  whereas  the  viscous  fingers 
were  damped  out  before  traveling  very  far  in  the  systems 
with  a  residual  water  saturation.  The  damping  mechanism  was 
attributed  to  the  movement  of  the  two  phases  in  the 
transverse  direction  to  flow. 

Hagoort  (1974)  also  studied  the  displacement  stability 
of  a  Buckley-Leverett  type  displacement.  Included  in  the 
study  were  the  effects  of  viscous  and  capillary  forces,  but 
gravity  was  neglected.  In  place  of  the  conventional 
end-point  mobility  ratio,  Hagoort  (1974)  defined  a  shock 
mobi 1 i ty  ratio  as : 


2.86 


It  was  found  that  the  displacement  is  unstable  provided  that 
the  shock  mobility  ratio  is  greater  than  unity,  subject  to 
the  constraint  that  the  wavelength  of  the  finger  is  smaller 
than  the  diameter  of  the  system.  It  was  further  concluded  by 
Hagoort  (1974)  that  a  Buckley-Leverett  type  displacement  is 
less  prone  to  become  unstable  than  a  Muskat  displacement. 

Several  researchers  in  this  laboratory  (Kloepfer,  1975; 
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Wiborg,  1976;  Baird,  1978)  conducted  tests  in  unconsolidated 
sandpacKs.  It  was  observed  that  a  dramatic  decrease  in 
recovery  occurred  at  high  rates.  This  decrease  was 
attributed  to  the  onset  of  instabilities.  However,  no 
physical  evidence  supporting  the  presence  of  viscous 
fingering  was  reported. 

As  recently  as  1979,  Peters  (1979)  investigated  the 
stability  problem  of  immiscible  displacements  in  systems 
with  and  without  a  residual  water  saturation.  Using  a  simple 
and  novel  technique,  he  was  able  to  capture  the  viscous 
fingers  and  obtain  permenent  records  via  the  use  of 
photography. 

Peters  (1979)  extended  the  theory  of  Chuoke,  et  al . 
(1959)  by  combining  all  of  the  pertinent  variables  of  the 
stability  problem  to  form  a  single  dimensionless  group.  He 
concluded  that  this  group  would  be  capable  of  predicting  the 
onset  of  instability.  Proceeding  in  a  manner  analogous  to 
Chuoke,  et  al .  (1959),  Peters  (1979)  showed  that  the 
stability  index  for  a  cylindrical  coordinate  system  may  be 
written  as: 


n 


2.87 


The  subscript  i  indicates  that  the  eigenvalues  can  take  on 
different  discrete  values.  Equation  2.87  is  simply  a 
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convenient  rearrangemnt  of  Equation  2.75.  As  before,  the 
necessary  and  sufficient  condition  for  instability  may  be 
written  as: 
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and  the  onset  of  instability  is  given  by: 


(ft  -  1 ) (V  -  V  )y  a2 
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2.89 


The  possible  discrete  values  of  6. a  are  given  by  the  zeros 
of  the  Bessel  function,  J '  (6 .  a)  .  The  first  non-zero  value 
of  6-a  is  1.8412  at  m=1,  and  consequently  the  onset  of 
instability  is  defined  as: 


(M  -  1  )(V  -  V  )y!  , 

r  c.  -w—  =  (1 .8412)^ 


Y*K 


2.90 
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Equation  2.90  may  be  written  as: 


(Mr  -  1)(V  -  Vc)ywD* 
C*^KWr 


13.56 


2.91 


where 


Q*  =  y*/y  2 . 92 

The  dimensionless  group,  given  by  Equation  2.91, 
defines  the  onset  of  instability  in  a  cylindrical  sandpack. 
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Provided  this  group  is  less  than  13.56,  the  displacement  is 
unconditionally  stable.  In  a  rectangular  coordinate  system, 
the  analogous  expression  to  Equation  2.91  is  given  by: 


-  1)(V  -  Vjy  L^L2 
c  w  x  y 


7T“ 


2.93 


It  should  be  noted  that  the  constant  of 
proportionality,  relating  the  bulk  interfacial  tension  and 
effective  interfacial  tension,  is  not  identical  to  that  of 
Chuoke,  et  al.  (1959).  With  this  definition,  Chuoke' s 
con  s  t  a  n  t  now  become  s : 


C  =  27t/3C* 


2.94 


From  the  photographic  records  of  the  viscous  fingers,  Peters 
(1979)  estimated  the  most  probable  wavelength.  Using 
Equation  2.85,  he  calculated  a  value  for  C  and,  hence,  C* 
(Equation  2.93)  and  y*  (Equation  2.92).  The  system 
parameters  obtained  by  Peters  (1979)  are  shown  in  Table  2.1. 

Table  2.1:  System  Parameters 

(Peters,  1979) 

Wettability  C  C*  y* 

dynes/cm 

Water-wet  190.45  306.25  7.442x103 

Oil-wet  25.40  5.45  1 . 324x 1 0 2 
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For  a  stability  number  greater  than  1000,  Peters  (1979) 
observed  that  the  recovery  was  insensitive  to  variations  of 
the  parameters  in  the  stability  number.  This  is  in  agreement 
with  the  observations  of  de  Haan  (1959).  At  a  stability 
number  greater  than  1000,  it  was  seen  by  Peters  (1979)  that 
numerous  fingers  with  a  small  wavelength  are  present.  The 
linear  theory  does  not  predict  this  second  breakpoint  in  the 
recovery  curve.  This  deficiency  may  be  attributed  to  the 
fact  that  linear  theory  is  inapplicable  for  disturbances 
with  amplitudes  much  larger  than  their  wavelength.  It  may 
also  be  because  the  current  theory  does  not  account  for 
viscous  dissipation  effects.  The  introduction  of  viscosity 
into  the  analysis  eliminates  the  tendency  for  the 
disturbances  of  small  wavelength  to  increase  without  bound 
(Bellman  and  Pennington,  1954). 

The  experimental  results  of  Peters  (1979)  agreed  well 
with  the  estimate  of  Chuoke,  et  al.  (1959).  However,  before 
the  stability  criterion  may  be  applied,  a  complete  set  of 
experiments  must  be  conducted  in  order  to  determine  the 
system  parameters.  As  a  consequence  of  the  extension  of 
Chuoke' s  theory,  Peters  (1979)  was  successful  in  predicting 
the  onset  of  instabilities.  His  experimental  estimates  of 
the  system  parameters  were  for  a  particular  fluid-rock 
system  and  a  complete  investigation  is  pending. 


3.  Experimental  Equipment,  Materials  and  Procedure 


3.1  Porous  Media 

3.1.1  Sand 

The  unconsolidated  sandpacks  used  throughout  this 
series  of  experiments  was  an  Ottawa  silica  sand  with  a  grain 
size  range  of  80-120  mesh  (Fisher  Scientific  S - 1 5 1 ) .  The 
sand  was  used  as  received  from  the  supplier  without  any 
screening . 

3.1.2  Coreholders 

The  sandpacks  were  housed  in  coreholders  constructed  of 
polyvinyl  chloride  (PVC)  for  low  pressure  usage  (less  than 
100  psi )  and  of  stainless  steel  for  high  pressure  studies. 
The  design  of  the  coreholders  was  the  same  as  those 
described  in  detail  by  Wiborg  (1976). 

Three  lengths  were  used  ranging  from  23  cm  to  110  cm 
with  a  choice  of  two  diameters  (2.4  cm  and  4.8  cm).  This  set 
of  coreholders  allowed  the  study  of  geometrically  dissimilar 
sandpacks . 

3.1.3  Packing  Procedure 

Several  techniques  were  used  in  packing  the  coreholders 
in  an  attempt  to  find  a  method  which  would  yield  sandpacks 
of  consistent  properties.  It  was  found  that  the  routine 
method  of  vibrating  a  dry  core  (Wiborg,  1976;  Baird,  1978) 
and  tamping  with  a  hammer  (Peters,  1979)  produced  sandpacks 
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with  properties  very  much  dependent  on  the  length  of  time  of 
vibration  or  tamping. 

Consequently,  the  coreholders  were  packed  with  wet  sand 
in  the  following  manner.  The  coreholder  was  partially  filled 
with  clear  distilled  water  followed  by  sand.  The  coreholder 
was  vibrated  throughout  the  filling  process  and  an 
approximately  constant  4  inch  head  of  water  was  Kept  in  the 
coreholder.  The  coreholder  was  then  vibrated  for  at  least  24 
hours.  After  the  vibration  period,  predried,  compressed  air 
at  room  temperature  was  passed  through  the  sandpack  for  48 
hours  to  ensure  the  complete  removal  of  the  packing  water. 
This  method  of  packing  is  a  modification  of  that  used  by 
Kloepfer  ( 1975) . 


3.2  Fluids 

Distilled  water  was  used  to  create  a  residual  water 
saturation,  and  distilled  water  doped  with  sodium 
fluorescein  was  used  as  the  displacing  fluid. 

Only  one  oil  was  used  during  the  tests.  This  was  a  Dow 
Corning  200  silicone-base  oil.  This  oil  was  used  because  its 
density  was  almost  equal  to  that  of  distilled  water.  The 
physical  properties  of  interest  for  these  two  fluids  are 
listed  in  Table  3.1  and  were  obtained  from  the  report  of 
Peters  ( 1 979 ) . 
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T ABLE :  3.1  Pertinent  Fluid  Properties  At  21.5*C 

Density  Viscosity 
(gm/cc)  (cp) 

Disti 1  led  Water  0.9910  1.028 

Dow  Corning  200  0.9667  105.4 

Interfacial  Tension:  24.3  dynes/cm 

Viscosity  Ratio:  102.5 

Water-oil  Density  Difference:  0.0243  gm/cc 


3.3  Pumping  System 

Two  dual  positive-displacement  RusKa  pumps  were  used  in 
this  investigation.  Both  pumps  were  capable  of  delivery 
rates  ranging  from  2.5  cc/hr  to  1120  cc/hr.  A  minor 
modification  to  the  external -drive  gearing  on  one  of  the 
pumps  was  done  in  order  to  attain  a  delivery  rate  of  0.5 
cc/hr.  The  modification  was  the  replacement  of  the  two  GB-60 
gears  by  a  GB-20  and  GB-100  gear.  Thus,  when  used  in  tandem, 
the  pumping  system  delivered  rates  from  0.5  cc/hr  to  2240 
cc/hr.  This  range  far  surpassed  that  of  other  investigators 
(Kloepfer,  1975;  Baird,  1978;  Peters,  1979) 

The  pumps,  charged  with  mercury,  displaced  the  fluid  of 
interest  from  a  stainless  steel  fluid  bomb  into  the 
coreholder.  All  lines  were  of  high-pressure  stainless  steel 
tubing.  The  complete  pumping  system  was  designed  to  a 
working  pressure  of  1800  psi  when  using  the  stainless  steel 
coreholders,  and  100  psi  for  the  PVC  coreholders.  Figure  3.1 
is  a  schematic  of  the  experimental  set-up. 
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SCHEMATIC  OF  DISPLACEMENT  APPRATUS 
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3.4  SandpacK  Property  Determination 

Prior  to  the  waterflood,  certain  properties  of  each 
sandpacK  were  measured.  These  measurements  were  performed  by 
routine  techniques  and  are  only  briefly  described. 

To  determine  the  pore  volume  (porosity),  the  coreholder 
was  first  evacuated  and  then  the  initial  fluid  (depending  on 
the  wettability  case  desired)  was  introduced  into  the 
sandpacK.  The  pore  volume  was  then  determined  by  both  a 
material  balance  calculation  and  by  weight.  Taking  into 
account  the  volume  of  fluid  contained  within  the  endcaps , 
agreement  was  within  2%. 

This  initial  fluid  was  then  flowed  through  the  sandpacK 
at  several  rates  in  order  to  estimate  the  absolute 
permeability  using  Darcy's  law.  In  the  case  of  water  only 
flowing,  a  mercury  manometer  was  used;  for  oil  only  flowing 
a  Heise  pressure  gauge  was  used. 

For  the  case  of  a  water-wet  system,  the  residual  water 
saturation  was  then  established  by  displacing  the  water  with 
oil.  The  displacement  was  conducted  at  a  high  rate,  with  the 
coreholder  in  the  vertical  position,  and  the  oil  was 
introduced  at  the  top. 

It  was  suggested  by  Saeedi  (1979)  that  the  rate  of 
injection  of  oil,  while  establishing  the  residual  water 
saturation,  should  be  higher  than  the  injection  rate  of  the 
subsequent  waterflood.  This  procedure  was  suggested  in  order 
to  ensure  that  the  residual  water  would  not  be  mobile.  With 
this  fact  in  mind,  injection  rates  were  set  so  that  the 


. 


. 


48 


maximum  pressure  drop  across  the  sandpack  did  not  exceed  100 
psi  ,  the  working  pressure  of  the  PVC  coreholder. 
Consequently,  this  restriction,  in  general,  was  sufficient 
to  satisfy  the  suggestion  of  Saeedi  (1979). 

Subsequently,  the  permeability  to  oil  at  a  residual 
water  condition  was  measured  in  the  same  manner  as  above.  It 
was  observed  that  no  water  or  sand  was  produced  during  this 
per i od . 


3.5  Displacement  Procedure 

After  the  sandpack  properties  were  measured,  the  flood 
was  conducted  in  the  following  manner.  The  coreholder  was 
placed  in  a  horizontal  position  inside  a  cabinet  maintained 
at  21.5’C.  The  inlet  end  was  connected  to  the  doped  water 
bomb  and  a  pressure  gauge.  A  graduated  cylinder  was  placed 
at  the  outlet  to  collect  the  effluent. 

After  an  appropriate  pumping  rate  was  selected,  the 
displacement  was  conducted  until  breakthrough  of  water. 
Breakthrough  was  that  event  when  the  first  drop  of  water  was 
produced.  The  use  of  the  coloured  water  aided  greatly  in 
determining  this  point  in  time.  At  this  time,  both  the  pump 
reading  and  volume  of  effluent  collected  were  recorded  in 
order  to  verify  the  absence  of  leaks  in  the  system, 
especially  in  the  tests  that  were  of  long  duration  (2-3 
weeks ) . 

Vertical  waterfloods  were  performed  in  the  same  manner 
as  the  horizontal  floods.  The  orientation  of  the  coreholder 
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was  with  the  inlet  end  down,  because  the  water  density  was 
slightly  higher  than  that  of  the  oil. 

At  the  end  of  the  test,  the  coreholder  was  disassembled 
completely,  cleaned,  and  repacked  with  fresh  sand  before 
commencing  another  test. 


3.6  Core  Extraction  Technique 

The  removal  of  the  sand  from  the  coreholder  was 
investigated  in  order  to  develop  a  technique  which  was 
quick,  simple  and  free  of  mess. 

The  method  normally  employed  in  this  laboratory  was  to 
use  a  long  auger  to  'drill7  out  the  core.  A  second  method 
was  to  use  hot  water  to  flush  out  the  sand  from  the 
coreholder.  Both  techniques  proved  to  be  time  consuming  and 
messy. 

The  following  procedure  was  developed.  After  the  outlet 
end-cap  was  removed,  the  coreholder  was  hung  in  a  vertical 
position  with  the  outlet  end  down.  The  Ruska  pumping  system, 
set  at  560  cc/hr,  was  used  to  inject  water  into  the 
coreholder.  This  had  the  effect  of  shifting  the  entire 
sandpack  several  centimeters.  After  the  initial  shift,  the 
waterline  was  replaced  by  a  90-psi  compressed  airline.  The 
compressed  air  ejected  the  entire  sandpack  from  the 
coreholder  into  a  bucket. 

This  technique  of  extracting  the  core  was  found  to  be 
both  rapid  (5-10  minutes)  and  relatively  free  of  mess. 
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4.  Presentation  and  Discussion  of  Results 
There  were  two  types  of  immiscible  displacement  tests 
carried  out  in  this  study.  The  first  type  was  conducted  in  a 
water-wet  porous  medium  with  a  residual  water  saturation. 

All  tests  of  this  type  were  horizontal  displacements. 

Secondly,  the  behaviour  of  an  oil-wet  porous  medium  was 
examined.  This  system  contained  no  residual  water  and  tests 
were  conducted  both  horizontally  and  vertically.  The 
fluid-rock  system  and  the  shape  of  the  coreholders  was 
identical  to  that  used  by  Peters  (1979)  in  order  to  make 
direct  use  of  his  data  base. 


4.1  Sandpack  Properties 

The  techniques  for  packing  each  coreholder  were 
outlined  in  Chapter  3.  It  was  observed  that  the  wet -packing 
method  produced  sandpacks  of  an  absolute  permeability  in  the 
range  of  14  to  19  Darcys.  However,  the  absolute  permeability 
of  the  dry,  vibrated  cores  was  generally  higher  and  quite 
unpredictable.  Wide  variations  in  permeability  were  seen, 
with  one  sandpack  having  a  permeability  of  32  Darcys. 

Similar  observations  were  made  for  the  permeability  to  oil 
at  residual  water  saturation. 

Both  the  porosity  and  residual  water  saturation 
appeared  to  be  insensitive  to  the  packing  procedure  used. 

The  porosity  of  the  sandpacks  was,  on  the  average,  of  the 
order  of  36%  and  the  residual  water  saturation  was 
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approximately  11%. 

Table  4.1  is  a  summary  of  the  sandpack  properties.  The 
order  of  appearance  of  the  data  does  not  reflect  the 
chronological  sequence  of  the  displacement  tests. 

The  lower  permeabilities  obtained  and  the  predictable 
results  of  the  wet  packing  method  may  be  attributed  to  the 
absence  of  static  electricity  and  the  stabilizing  effect  of 
water.  The  inherent  sorting  phenomenon  associated  with 
vibrating  is  consequently  suppressed.  From  these 
observations,  it  may  be  concluded  that  the  wet -packing 
method  is  superior  to  other  techniques. 

An  experiment  was  designed  in  order  to  verify  the 
complete  removal  of  water  from  the  sandpack  during  the 
packing  procedure.  The  coreholder  was  prepared  dry,  then 
saturated  and  dryed  with  compressed  air.  It  was  observed 
that,  based  on  measurements  of  the  weight  of  the  coreholder, 
all  the  water  was  removed  by  the  drying  procedure. 

Although  no  prescreening  of  the  sand  was  attempted,  it 
is  suggested  that  a  complete  set  of  tests  be  arranged  with 
sand  from  the  same  lot.  It  was  observed  that  a  change  in 
lots  corresponded  to  the  dramatic  change  in  permeability 
seen  in  Table  4.1  (Runs  18,  34,  and  35). 


4.2  Wettability  of  the  Porous  Media 

When  clean  and  uncontaminated,  Ottawa  silica  sand  is 
preferentially  water-wet.  However,  the  wettability  of  porous 
media  depends  on  both  the  oil  composition  and  rock  type. 


' 


. 


52 


TABLE:  4.1 

SUMMARY  OF  SANDPACK  PROPERTIES 
WATER-WET  SYSTEM 


RUN 

LENGTH 

DIAMETER 

POROSITY 

ABS.  PERM. 

PERM.  @SWI 

SWI 

(CM) 

(CM) 

(%) 

(DARCY) 

(DARCY) 

(%) 

ID 

113.1 

2.390 

36.36 

21.04 

17.25 

10.4 

2D 

113.1 

2.390 

37.74 

19.38 

17.07 

10.4 

3D 

113.1 

2.390 

37.  19 

19.72 

14.61 

11.5 

4W 

113.1 

2.390 

37.92 

18.34 

15.30 

12.6 

5W 

26.9 

2.390 

38.01 

18.26 

16.59 

10.8 

6W 

22.8 

4.859 

34.10 

14.23 

12.87 

10.3 

7W 

22.6 

4.815 

37.28 

17.91 

15.28 

10.2 

8W 

22.8 

4.859 

33.66 

16.60 

11.67 

9.4 

9W 

53.2 

4.850 

32.32 

15.38 

14.90 

10.7 

10W 

53.2 

4.850 

36.98 

18.21 

14.93 

10.7 

1 1 W 

53.2 

4.850 

35.83 

18.37 

14.51 

13.2 

1 2D 

53.2 

4.850 

37.56 

20.83 

15.88 

12.4 

1  3D 

110.4 

4.864 

33.65 

19.  17 

16.98 

10.7 

14W 

110.4 

4.864 

36.22 

18.35 

14.48 

11.9 

OIL- 

WET  SYSTEM 

RUN 

LENGTH 

DIAMETER 

POROSITY 

ABS.  PERM. 

(CM) 

(CM) 

(%) 

(DARCY) 

15D 

113.1 

2.390 

37.29 

18.83 

16D 

113.1 

2.390 

37.98 

19.48 

17D 

113.1 

2.390 

37.09 

19.90 

18W 

26.9 

2.390 

38.77 

10.81 

19D 

22.6 

4.815 

35.50 

17.  15 

20D 

22.8 

4.859 

35.97 

16.93 

2 1 W 

53.2 

4.850 

35.25 

17.70 

22W 

53.2 

4.850 

35.60 

17.45 

23W 

53.3 

4.811 

35.80 

18.02 

24W 

53.3 

4.811 

36.69 

18.57 

25W 

53.3 

4.81 1 

38.58 

21.98 

26D 

110.4 

4.864 

38.41 

22.83 

27D 

110.2 

4.828 

36.67 

19.78 

28D 

110.2 

4.828 

37.  13 

18.88 

29D 

110.2 

4.828 

39.20 

31.70 

30D 

110.2 

4.828 

36.79 

18.91 

3 1 D 

110.4 

4.864 

36.63 

19.51 

Continued  Next  Page 
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TABLE:  4. 1  CONTINUED 
SUMMARY  OF  SANDPACK  PROPERTIES 

OIL-WET  SYSTEM- -VERT ICAL  FLOODS 


RUN 

LENGTH 

(CM) 

DIAMETER 

(CM) 

POROSITY 

(%) 

ABS.  PERM 
(DARCY) 

32W 

22.8 

4.859 

36.97 

18.56 

33W 

22.6 

4.815 

37.20 

17.84 

34W 

22.6 

4.815 

38.  12 

10.75 

35W 

22.8 

4.859 

35.74 

10.76 

D 

W 


Dry  Packed 
Wet  Packed 
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The  effect  of  Dow  Corning  200  fluid  on  the  wettability 
of  Ottawa  silica  sand  was  investigated  by  Baird  (1978).  For 
a  water-wet  system,  his  data  suggested  that  Dow  Corning 
fluid  tends  to  shift  the  wettability  of  the  Ottawa  silica 
sand  toward  oil -wet.  Thus,  the  Dow  Corning  fluid  tends 
preferentially  to  wet  the  sand. 

Based  on  Baird's  (1978)  findings,  the  porous  media  in 
this  study  were  considered  to  be  wet  by  the  fluid  of  first 
contact.  This  is  provided  that  the  oil  was  not  permitted  to 
remain  in  contact  with  the  sand  for  a  long  period  of  time. 
Consequently,  for  a  displacement  test  where  water  was  the 
first  fluid  to  contact  the  sand,  the  system  was  considered 
to  be  water-wet.  Furthermore,  the  porous  media  were 
considered  to  be  oil-wet  when  Dow  Corning  200  was  the 
initial  saturating  fluid.  It  was  beyond  the  scope  of  this 
study  to  establish  a  residual  water  saturation  in  the  latter 
system. 


4.3  Displacement  Results 

Immiscible  displacement  tests  were  carried  out  for  both 
wettability  types  over  a  wide  range  of  rates  varying  from 
0.5  cc/hr  to  480  cc/hr .  The  coreholders  were  of  three 
lengths  and  two  diameters.  Due  to  the  small  pore  volumes 
associated  with  the  2.4  cm  coreholders,  the  flow  rates  were 
restricted  to  a  maximum  of  140  cc/hr  for  the  longer  core  and 
to  40  cc/hr  for  the  shorter- length  sandpacKs . 
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Table  4.2  is  a  summary  of  the  results  of  the 
displacement  tests  carried  out  in  this  study.  Included  in 
the  table  are  the  results  of  the  calculations  of  the  various 
similarity  groups  pertinent  to  this  study. 


4.4  Scaling  in  Porous  Media--Theoretical  Results 

In  a  scaled  model  study,  it  is  not  necessary  to  demand 
that  the  geometric  factor  be  identical  in  both  model  and 
prototype,  provided  that  the  displacement  is  stable. 

However,  the  model  and  prototype  should  have  the  same  shape 
because  the  derivation  of  the  scaling  groups  was  based  on  a 
particular  coordinate  system. 

The  meaning  of  the  terms  'geometric  factor'  and 
'geometric  similarity'  should  not  be  confused.  Geometric 
factor  refers  to  the  ratio  of  the  overall  dimensions  (L/D) 
of  the  porous  medium  whereas  geometric  similarity  refers  to 
the  scaling  of  the  characteristic  lengths  of  the  porous 
med  i  urn . 

The  mul ti -dimensional  modified  inspect ional  analysis 
results  in  the  emergence  of  two  capillary  numbers  in  the  set 
of  similarity  groups.  This  important  result  does  not  arise 
if  the  displacement  problem  is  assumed  to  be  linear.  The 
linear  capillary  number  dictates  the  stabilization  of  the 
front  in  the  direction  of  bulK  flow.  The  radial  capillary 
number,  which  is  important  only  if  an  appreciable  saturation 
gradient  exists  in  the  radial  di rect ion--such  as  in  the  case 
of  viscous  f i nger i ng- -control s  the  spreading  or  growth  of 
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the  'front'  in  the  radial  direction. 

Although  not  presented,  the  similarity  groups  which 
arise  from  a  three-dimensional  analysis  in  a  Cartesian 
coordinate  system  include  three  capillary  numbers.  These  may 
be  interpreted  in  the  same  way  as  the  linear  and  radial 
capillary  numbers. 


4.5  Scaling  in  Porous  Medi a--Exper imental  Results 

The  breakthrough  recovery  was  correlated  with  both  the 
conventional  similarity  groups  which  arise  from  an 
inspectional  analysis  using  the  Leverett  J-function 
(Appendix  A. 2. 2),  and  those  derived  in  this  study  (Appendix 
A.2.3).  For  each  wettability  type  considered,  it  was  assumed 
that  the  mobility  ratio,  relative  permeability  curves,  and 
capillary  pressure  curves  remained  constant.  Furthermore,  it 
was  assumed  that  gravity  effects  were  negligible.  The 

similarity  groups  of  interest  were: 

1  /  2 

1)  y(K<J>)  '  /VuwL  Linear  Leverett  number  Lt£ 

and  its  counterpart, 

crK  /v u  L  Linear  capillary  number  N  0  ; 

w  r  w  c * 


2) 


y(K<j))^ ^/Vy  D  Radial  Leverett  Number  Ltr 

and  its  counterpart, 

aK  /Vu  D  Radial  capillary  number  N  ; 

wr  w  cr 


3)  L/D 


Geometric  factor. 
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4.5.1  Water-wet  System 

In  this  type  of  system,  the  breakthrough  recovery,  when 
correlated  with  the  inverse  linear  Leverett  number, 
exhibited  the  behaviour  shown  in  Figure  4.1.  It  was 
observed,  as  predicted  by  theory,  that  the  breakthrough 
recovery  is  insensitive  to  variations  in  the  similarity 
group,  provided  that  the  displacement  is  stable  and 
stabilized.  This  region  is  represented  by  the  horizontal 
line  (and  the  dashed  portion  is  an  extension  according  to 
theory)  at  a  recovery  of  41%  IOIP.  Furthermore,  regardless 
of  the  sandpack  dimensions,  all  data  correlated  well.  Thus 
the  recovery  is  insensitive  to  the  geometric  factor  in  this 
region . 

However,  a  deviation  from  this  plateau  is  very  apparent 
at  higher  values  of  the  inverse  linear  Leverett  number.  This 
may  be  attributed  to  the  viscous  forces  beginning  to 
dominate  the  recovery  process.  The  Points  A,  B,  C  and  D 
correspond  to  the  onset  of  instability  at  the  floodfront  for 
the  various  geometric  factors  considered  in  this  study. 

These  points  correspond  to  a  stability  number  of  13.56  and 
were  determined  experimentally. 

For  unstable  displacements,  a  family  of  curves  results. 
As  shown  in  Figure  4.1,  each  member  of  the  family 
corresponds  to  a  particular  L/D  ratio.  Thus,  it  appears  that 
the  geometric  factor  cannot  be  omitted  as  a  scaling  group 
when  viscous  fingering  is  present.  Consequently  in  this 


. 


FIG.  4.1:  RECOVERY  CORRELATED  BY 
THE  INVERSE  LINEAR  LEVERETT  NUMBER 
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intermediate  zone,  where  the  viscous  forces  begin  to 
dominate,  the  linear  flow  assumption  is  no  longer  valid. 
Hence,  scaling  groups  based  on  the  assumption  of  a  stable 
displacement  may  not  be  valid  if  instabilities  are  present, 
although  they  may  still  be  useful  if  correctly  interpreted. 

Based  on  the  results  of  Baird  (1978),  a  second  region 
(which  may  be  considered  as  a  pseudo-stable  displacement 
region)  of  constant  recovery  is  postulated.  This  region  is 
represented  by  the  lower  dashed  curve  in  Figure  4.1.  Points 
E,  F,  G  and  H  correspond  to  where  the  viscous  forces 
completely  dominate  the  recovey  process,  that  is,  for  a 
stability  number  of  900.  These  points  were  calculated  using 
estimated  average  sandpacK  and  fluid  properties.  A  typical 
calculation  is  presented  in  Appendix  C.  It  may  be  speculated 
that  beyond  these  points,  the  recovery  is  again  insensitive 
to  the  magnitude  of  the  parameters  in  the  linear  Leverett 
number.  This  speculation  is  based  on  the  results  of  the 
oil-wet  tests,  which  are  presented  in  a  later  section.  It 
should  be  noted  that  there  is  little  experimental  basis  for 
the  level  at  which  the  lower  dashed  portion  of  the  recovery 
curve  has  been  plotted.  As  will  be  seen  in  a  later  section, 
the  level  of  the  recovery  is  partially  justified  by  the 
extrapolation  of  the  recovery  curve  to  a  stability  number  of 
900. 

The  pseudo- stable  region  could  not  be  obtained 
experimentally  due  to  the  high  displacement  rates  required, 
approximately  3200  cc/hr ,  and  the  small  pore  volumes  of  the 


. 
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coreholders.  The  largest  sandpack,  with  an  estimated 
recovery  of  17%  I 0 I P ,  would  break  through  in  less  than 
1  minute  at  this  rate.  Consequently,  data  from  this  type  of 
test  were  considered  invalid. 

The  correlation  of  the  breakthrough  recovery  with  the 
radial  Leveret t  number  (shown  in  Figure  4.2)  is  acceptable, 
provided  the  sandpack  diameters  are  approximately  equal. 
Although  it  appears  that  the  data  correlate  well,  regardless 
of  geometric  factor,  it  should  be  noticed  that  Points  A  (for 
the  4.8  cm  diameter  sandpack)  and  B  (2.4  cm  diameter 
sandpack)  both  correspond  to  the  onset  of  frontal 
perturbations.  This  'good  correlation'  can  be  explained  by 
the  relatively  small  difference  in  magnitude  of  the  sandpack 
diameters,  and  by  the  wide  scatter  of  the  data  in  the 
perturbed  region. 

Peters  (1979)  suggested  that  the  radial  Leverett  number 
is  the  better  scaling  parameter  in  systems  which  contain  a 
residual  water  saturation.  This  conclusion,  however,  was 
based  on  sandpacks  of  similar  diameter.  Despite  the 
difficulty  of  the  interpretation  for  the  two  reasons 
mentioned  above,  the  data  do  not  correlate  well  for 
different  diameters. 

From  the  above  discussion,  it  can  be  concluded  that 
neither  the  linear  Leverett  number  nor  the  radial  Leverett 
number  is  a  good  correlating  group  beyond  the  point  where 
viscous  forces  begin  to  dominate  the  process.  Thus,  a  model 
may  be  compared  to  a  prototype  provided  the  displacement  is 
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FIG.  4.2:  RECOVERY  CORRELATED  BY 
THE  INVERSE  RADIAL  LEVERETT  NUMBER 


63 


w00Q  £>cn 
cc.^%  +<J) 


lu 

2:  ^  _  CO 

a:  •©  ❖cc 

LU 

CD  I— 


LD 

OJ 

o 


o 


cocoox"  a: 

CM  LO  — •  Q- 

>- 


o 

I — 

o 

z 

LU 


o 

cc 


CL  CD 
I—  o 
a:  cj 
a  co 

O  1 — 1 

O  > 


o 

o 


T 

o 

cn 


T 

O 

00 


T 

o 

r*'* 


LUi 


T 

O  O 

CD  LO 


O 

=3* 


o 

o 


-  -  CM 


--u> 


CD 


a 


o 

CD 


•e- 


CJ 

I 

0 


02 

UJ 

00 


--  CJ 


UJ 

02 

LU 


6  ^ 


in 


CJ 


CO 

CE 

02 

a*  LU 

0 

1  02 

LU 


m 


CJ 


"I - 1 - T 

0000 


in 

1 

CD 


CO  C\J  ri 

t d i o i z.)  Ay3A033y  HongyHi>iH3y8 


64 


stable.  Moreover,  scaling  the  geometric  factor  is  not 
important  for  such  displacements. 

Figures  4.3  and  4.4  are  similar  to  the  previous  two 
figures,  except  that  the  correlation  of  the  breakthrough 
recovery  is  made  with  the  linear  capillary  and  radial 
capillary  numbers,  respectively.  It  was  assumed  that  the 
permeability  to  water  at  residual  oil  saturation  was  equal 
to  the  absolute  permeability.  The  effect  of  this  assumption 
is  to  shift  the  curve  horizontally  to  the  left.  A  test  was 
undertaken  in  which  the  sandpack  was  flooded  out  to  a 
water-oil  ratio  of  49.  The  permeability  to  water  at  residual 
oil  saturation  was  estimated  to  be  6  Darcys.  This  result  is 
in  agreement  with  the  experiments  of  Saeedi  (1979).  It  is 
possible  to  assume  that  the  permeability  to  water  at 
residual  oil  saturation  is  some  fixed  fraction  of  the 
absolute  permeability.  However,  this  entails  the  implicit 
assumption  that  the  ratio  is  fixed  throughout  the  set  of 
tests.  No  experimental  evidence  is  available  to  support  this 
assumption.  For  the  sake  of  consistency,  the  absolute 
permeability  was  used  in  all  calculations. 

Both  Figures  4.3  and  4.4  were  constructed  assuming  the 
capillary  pressure  normalizing  parameter  to  be  of  the  order 
of  2 . 027x 1 04  dynes/cm2  (Saeedi,  1979).  The  validity  of  this 
assumption  was  verified  by  equating  the  Leveret t  number  and 
capillary  number,  that  is: 

v  /Kd>  _  C7  K 

V  -  v^u 


4.  1 


FIG.  4.3:  RECOVERY  CORRELATED  BY 
THE  INVERSE  LINEAR  CAPILLARY  NUMBER 
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or 


a  =  y/<J>'/K  4.2 

Using  average  properties,  this  estimation  showed  that  the 
capillary  pressure  normalizing  parameter  is  of  the  order  of 
3.5x104  dynes/cm2. 

Thus,  the  correlation  with  the  capillary  number  shows 
that  all  the  displacement  tests  carried  out  in  this  study 
were  unstabilized  as  defined  in  Chapter  1,  that  is,  the 
inverse  linear  capillary  number  was  much  less  than  100.  This 
conclusion  could  not  be  drawn  by  past  researchers  for  the 
following  reasons. 

1)  No  estimate  of  the  magnitude  of  the  Leverett 
number  delimiting  a  stabilized  displacement  was 
avai 1  able . 

2)  The  plateau  observed  in  the  Leverett  correlation 
may  be  mistaken  for  a  stabilized  displacement 
wherein  the  results  may  be  masked  by 
'end-effects'  (Kyte  and  Rapoport,  1958).  The 
observed  plateau  is  due  to  an  excessive  build  up 
of  water  at  the  outlet  end.  Consequently,  the 
observed  recovery  is  too  high. 

3)  Figure  2.8  shows  a  typical  recovery  response 
curve.  It  can  be  seen  that  at  a  high  mobility 
ratio  and  large  values  of  the  capillary  number, 
the  recovery  is  affected  minimally  by  the 
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magnitude  of  the  parameters  in  the  Leverett 
number.  Thus,  the  recovery  is  only  a  weak 
function  of  the  Leverett  number,  provided  the 
mobi 1 i ty  ratio  is  high. 

Figures  4.1  (a  correlation  of  the  recovery  with  the  inverse 
Leverett  number)  and  4.3  (a  correlation  of  the  recovery  with 
the  inverse  capillary  number)  are  virtually  identical.  As  a 
consequence,  at  least  for  this  experimental  system,  the 
Leverett  number  and  capillary  number  have  the  same  physical 
significance.  This  suggests  that  the  magnitude  of  the  linear 
Leverett  number  should  be  of  the  order  of  0.01  if  the 
displacement  is  to  be  stabilized.  In  light  of  point  3  above, 
'stabilization'  may  be  defined  as  that  point  where  the 
recovery  response  curve  is  within  5%  of  the  steady-state 
recovery  as  predicted  by  Buck  ley- Leveret t  theory.  However, 
it  may  not  be  convenient  to  use  this  definition  due  to 
difficulties  in  determining  the  Buck  ley- Leveret t  recovery. 

4.5.2  Oil -wet  System- -Hor izont a  1  Displacements 

The  breakthrough  recovery,  for  the  oil -wet  system,  was 
correlated  against  the  same  scaling  groups  as  for  the 
water-wet  system.  Figure  4.5  shows  the  recovery  correlated 
with  the  inverse  linear  Leverett  number.  The  dashed  curves 
in  this  figure  indicate  the  'expected'  response  for  this 
type  of  system.  Again,  as  for  the  water-wet  case,  the  only 
experimental  basis  for  the  level  of  the  recovery  curve  is 
from  the  extrapolation  of  the  recovery  curve  to  a  stability 
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FIG.  4.5:  RECOVERY  CORRELATED  BY 
THE  INVERSE  LINEAR  LEVERETT  NUMBER 
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number  of  13.56.  Points  A,  B,  and  C  define  the  onset  of 
frontal  perturbations  while  Points  E,  F,  and  G  indicate  the 
complete  dominance  of  the  viscous  forces.  The  'expected' 
response  may  be  justified  by  the  same  argument  presented  for 
the  water-wet  system. 

In  an  attempt  to  explain  the  deviation  from  the 
'expected'  response,  the  sandpacKs  were  extracted  from  the 
coreholders  and  a  post-mortem  was  performed.  It  was  observed 
that  the  waterflood  suffered  from  gravity  underride  by  the 
injected  water  at  low  flow  rates.  It  should  be  noticed  that 
the  effects  of  gravity  are  more  pronounced  in  the  longer 
sandpacKs.  Furthermore,  the  experimental  results  are  masked 
by  the  compounded  effects  of  frontal  perturbations  and  a 
non-stabi 1 ized  displacement.  These  compounded  effects 
confuse  the  interpretation  of  the  data.  Any  interpretat ion 
must  be  carefully  undertaken. 

In  view  of  the  fact  that  a  water  tongue  was  observed, 
the  waterflood  at  low  displacement  rates  was  considered  to 
behave  according  to  the  theory  of  Dietz  (1953).  Based  on 
this  theory,  Hawthorne  (1960)  derived  the  following  equation 
to  calculate  the  slope  of  the  fluid-fluid  interface,  viz: 


tan  to 


Vcosa  (^w/^wr  "  uc/^0r^  _ 

-  P  }  +  Vsina  (yw/Kwr  -  M0/Kor) 


4.3 


where  the  angles  a  and  to  are  the  dip  of  the  formation  and 
interface  with  the  horizontal,  respectively.  For  a 
horizontal  displacement  and  assuming  that  the  end-point 


. 


permeabilities  are  equal  to  the  absolute  permeability, 
Equation  4.3  may  be  written  as: 
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4.4 


where 


u 


w 


0 


4.5 


Equation  4.4  may  be  rewritten  as: 


tan  03 


Ay 


4.6 


Using  Equation  4.6  and  the  data  from  Run  16,  the  angle  of 
the  interface  was  estimated  as  60°.  The  post-mortem  of  this 
sandpacK  enabled  the  angle  to  be  estimated  as  approximately 
55°,  in  fair  agreement  with  the  calculated  angle. 

The  results  of  the  water-wet  tests  showed  that  no 
correlation  of  all  the  data  was  possible  in  the  intermediate 
zone  (where  the  viscous  forces  begin  to  dominate  the 
displacement  process).  It  is  also  expected  that  this  should 
be  true  in  the  oil -wet  system.  However,  due  to  the  effects 
of  gravity,  the  data  in  this  intermediate  zone  tend  to  fall 
along  the  same  curve.  This  implys  that  a  correlation  is 
possible.  However,  the  conclusion  that  the  linear  capillary 
number  is  a  good  correlating  parameter  in  an  oil -wet  system 
is  incorrect. 
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The  existence  of  a  second  stabilized  zone  is  clearly 
observed  in  Figure  4.5.  In  this  region,  all  data  correlated 
well,  suggesting  that  the  geometric  factor  does  not  affect 
the  recovery. 

Figure  4.6  is  similar  to  Figure  4.2  and  the  same 
comments  apply.  Again,  the  effects  of  gravity  are  large  for 
slow  velocities  and  long  cores.  The  dashed  portion  of  the 
curve  indicates  the  'expected'  recovery,  if  these  effects 
are  negligible. 

Using  the  estimate  of  the  capillary  pressure 
normalizing  parameter,  Figures  4.7  and  4.8  were  constructed. 
It  should  be  noted  that  the  capillary  pressure  normalizing 
parameter  may  be  very  different  in  an  oil-wet  system,  as 
compared  to  a  water-wet  one.  Typically,  the  capillary 
pressure  normalizing  parameter  is  one  to  two  orders  of 
magnitude  smaller  in  an  oil -wet  system. 

The  breakthrough  recovery,  when  plotted  against  the 
inverse  gravity  number  (Figure  4.9)  is  remarkably  similar  in 
character  to  the  radial  scaling  groups.  It  can  be  seen  that, 
for  small  values  of  the  inverse  gravity  number,  the  recovery 
is  dramatically  reduced.  No  interpretation  is  possible  from 
the  data  at  higher  values  of  the  inverse  gravity  number  due 
to  the  effects  of  instabilities.  However,  it  may  be  stated 
that  the  effects  of  gravity  become  negligible  at  high 
displacement  rates. 
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FIG.  4.6:  RECOVERY  CORRELATED  BY 
THE  INVERSE  RADIAL  LEVERETT  NUMBER 
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FIG.  4.7:  RECOVERY  CORRELATED  BY 
THE  INVERSE  LINEAR  CAPILLARY  NUMBER 
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FIG.  4.8:  RECOVERY  CORRELRTED  BY 
THE  INVERSE  RRDIRL  CRPILLRRY  NUMBER 
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4.5.3  Oil-wet  System--Ver t ical  Displacements 

In  an  attempt  to  obtain  the  'expected'  recovery  shown 
in  Figure  4.5,  several  vertical  displacements  were 
conducted.  In  a  vertical  flood,  the  effects  of  gravity 
underride  are  negated  and  consequently,  the  displacement 
should  proceed  as  a  plane  horizontal  front. 

It  was  observed  that  all  the  vertical  tests  resulted  in 
recoveries  well  below  the  'expected'  recovery.  In  order  to 
gain  insight  into  the  apparent  failure  of  the  tests,  an 
examination  was  conducted  on  the  sandpacKs.  The  following 
observations  were  made. 

Run  32 

This  test  was  conducted  in  the  short  (4.8  cm  diameter,  23  cm 
length)  coreholder.  It  was  seen  that  the  areal  sweep 
efficiency  was  high  up  to  the  last  7  cm  of  the  sandpack. 
Between  this  point  and  the  outlet,  the  areal  sweep 
efficiency  decreased  rapidly.  The  sandpack  was  sliced  in 
1  cm  sections.  This  revealed  an  almost  circular  water 
channel  centred  in  the  sandpack,  the  diameter  of  the  circle 
decreasing  toward  the  outlet.  This  peculiar  behaviour,  near 
the  outlet,  was  attributed  to  the  design  of  the  end-cap.  It 
was  thought  that  the  end-cap's  effect  was  similar  to  that  of 
a  point  sink,  thus  distorting  the  flow  channel. 

Run  33 


This  test  was  conducted  in  the  same  manner  as  the  previous 


. 
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test,  except  for  a  modification  to  the  design  of  the  end-cap 
in  an  attempt  to  'remove'  the  sink  effect.  The  result  was  to 
increase  the  recovery  by  approximately  10%.  However,  this 
test  again  fell  short  of  the  'expected'  85%  recovery.  The 
post-mortem  showed  that  the  horizontal  interface  became 
distorted  within  4  cm  of  the  outlet.  The  effect  of  the 
end-cap  was  again  influencing  the  flow  path  and  a  sectioning 
of  the  sandpack  revealed  an  elliptically  shaped  flow  channel 
open  to  water.  Thus,  the  modification  was  only  partially 
successful . 

Run  34 

This  third  test  was  conducted  with  a  lower  permeability,  all 
other  parameters  being  the  same.  As  in  the  previous  tests, 
this  displacement  was  conducted  at  a  velocity  below  the 
critical  in  order  to  avoid  viscous  fingering.  The 
post-mortem  of  this  sandpack  showed  that  the  flow  path 
followed  by  the  water  was  cylindrical  in  nature,  that  is, 
the  flow  tended  to  follow  the  perimeter  of  the  coreholder. 

The  results  of  these  three  tests  suggest  that  the 
behaviour  of  the  waterflood  is  extremely  sensitive: 

1)  to  any  heterogeneities  in  the  porous  medium;  and 

2)  to  the  'boundary  conditions'  of  the  system 
(end-caps,  walls  of  the  coreholder) 

if  the  superficial  velocity  is  small. 
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Run  35 

This  test  was  conducted  in  the  regime  of  unstable 
displacements.  In  addition  to  the  effects  mentioned  above, 
it  was  observed  that  the  flow  paths  were  irregular. 


4.6  Stability  of  Displacement  Tests 

Displacement  tests,  for  both  types  of  wettability,  were 
conducted  using  sandpacKs  of  various  dimensions.  It  was  seen 
in  these  series  of  tests  that  a  displacement  may  fall  into 
one  of  the  following  regimes: 

1)  stable  displacement:  no  viscous  fingers  exists 
and  there  is  a  balance  between  viscous  and 
capi 1 1  ary  forces ; 

2)  unstable  displacement  (transition  zone):  viscous 
forces  begin  to  dominate  the  displacement  and 
the  displacement  is  characterized  by  a  sharp 
decrease  in  recovery; 

3)  unstable  displacement  (pseudo-stable):  the 
displacement  is  completely  dominated  by  viscous 
forces  and  the  recovery  is  insensitive  to 
instabi 1 i ties . 

4.6.1  Water-wet  System 

Figure  4.10  is  a  graphical  summary  of  the  stability  of 
the  water-wet  displacement  tests  performed  on  sandpacKs 
having  several  geometric  factors.  It  should  be  remembered 
that  the  stability  number  does  not  attempt  to  'scale'  the 
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displacement  performance.  Rather,  it  delimits  the  stability 
of  the  displacement  process. 

An  objective  of  the  study  was  to  define  better  the 
onset  of  instability.  To  this  end,  displacement  tests  were 
carried  out  in  the  vicinity  of  a  stability  number  of  13.56. 
As  predicted  by  theory,  the  data  from  this  study  supports 
the  findings  of  Peters  (1979).  It  was  observed  that  the 
recovery  for  this  particular  fluid-rock  system  was 
approximately  41%  IOIP,  provided  that  the  displacement  was 
stable.  Beyond  the  onset  of  instability  at  the  floodfront,  a 
sharp  decrease  was  observed,  corresponding  to  the  appearence 
of  viscous  fingers. 

Due  to  the  experimental  limitations  (high  displacement 
rates  required  and  small  pore  volumes  of  the  sandpacks),  it 
was  not  feasible  to  conduct  tests  beyond  a  stability  number 
of  400.  Thus,  it  was  not  possible  to  define  the  second 
breakpoint  in  the  recovery  curve.  In  order  to  define  the 
pseudo-stable  region,  a  more  viscous  oil  is  required  because 
an  increase  in  the  mobility  ratio  has  the  effect  of  shifting 
the  curve  down  and  to  the  right. 

The  data  base  of  Peters  (1979)  was  extended  by 
approximately  one  order  of  magnitude,  to  small  values  of  the 
stability  number.  The  displacement  tests  showed  that  the 
recovery  is  independent  of  a  stability  number,  provided  it 
is  less  than  its  critical  value,  13.56.  This  portion  of  the 
curve  corresponds  to  the  nearly  stabilized  displacements 
discussed  in  Section  4.5.1. 


- 

' 


82 


4.6.2  Oil -wet  System 

Figure  4.11  shows  the  scaled  onset  of  instability  for 
the  oil -wet  system  studied.  It  was  possible  to  extend  the 
data  base  of  Peters  (1979)  and  obtain  data  in  all  three 
regions,  as  defined  in  the  introductory  comments  of  this 
section.  The  dashed  portion  of  the  curve  indicates  the 
'expected'  recovery  as  discussed  in  an  earlier  section.  The 
vertical  displacement  results  are  also  indicated. 

According  to  the  observations  made  in  the  water-wet 
system,  all  data  should  plot  as  one  curve.  However,  the 
family  of  curves  indicated  in  Figure  4.11  are  a  consequence 
of  the  effects  of  gravity  underride.  It  should  be  noted  that 
the  gravity  effects  are  stronger  in  the  longer  length 
sandpacks.  The  effects  of  gravity  are  largest  in  the  stable 
displacement  region  and  lessen  with  increasing  stability 
number  (Note:  increasing  the  stability  number  corresponds  to 
increasing  the  superficial  velocity). 

At  a  stability  number  of  900,  all  data  converge.  It  was 
mentioned  in  Chapter  2  that,  as  yet,  there  is  no  theory 
available  to  predict  the  onset  of  this  pseudo-stable 
displacement  region. 

Beyond  a  stability  number  of  900,  the  recovery  is 
independent  of  the  stability  number.  Peters  (1979)  suggested 
that,  in  this  region,  the  recovery  should  decline 
moderately.  However,  the  additional  data  from  this  study 
indicates  that  the  recovery  is  independent  of  the  stability 
number  in  this  region. 
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5.  Summary  and  Conclusions 

The  conclusions  reached  from  this  study  may  apply  only 
to  the  particular  geometry  and  fluid-rock  system  used  here. 

A  broad  application  of  the  results  should  be  avoided  until  a 
more  thorough  investigation  is  undertaken. 

A  mul ti -dimensional  modified  inspectional  analysis  was 
undertaken  in  order  to  determine  a  set  of  similarity  groups 
pertinent  to  the  immiscible  fluid  displacement  problem.  The 
theoretical  analysis  led  to  the  following  conclusions. 

1.  The  set  of  similarity  groups  which  arise  from  a 

mul ti -dimensional  modified  inspectional  analysis  is: 
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for  a  cylindrical  coordinate  system.  Similarly,  in  a 
Cartesian  coordinate  system,  the  following  similarity 
groups  arise: 


KwrAp9 
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or  Mw 


2.  The  analysis  suggests  that  it  may  not  be  necessary  to 
demand  that  the  geometric  factor  be  identical  in  both 
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the  model  and  the  prototype,  provided  that  the 
displacement  is  stable. 

As  a  prelude  to  the  experimental  verification  of  the 
theoretical  findings  above,  the  core  packing  procedure  was 
investigated.  It  was  concluded  that  the  wet  packing 
technique  is  effective  in  preparing  sandpacks  of  predictable 
and  consistent  properties.  It  should  also  be  mentioned  that, 
as  gained  through  experience,  a  complete  series  of  tests 
should  be  arranged  with  sand  from  a  particular  lot  shipment. 

The  experimental  part  of  this  study,  with  respect  to 
the  scaling  of  the  displacement  problem,  led  to  the 
conclusions  below. 

1.  Neither  the  linear  nor  the  radial  scaling  groups  are 
good  correlating  parameters  beyond  the  point  where 
viscous  forces  begin  to  dominate  the  displacement 
process.  This  conclusion  is  valid  for  both  wettability 
types  studied.  Thus,  a  model  may  be  compared  to  a 
prototype,  provided  the  displacement  is  stable. 

Moreover,  scaling  the  geometric  factor  is  not  important 
for  such  displacements. 

2.  Scaling  groups,  based  on  the  assumption  of  a  stable 
displacement,  may  not  be  valid  if  instabilities  are 
present.  However,  they  do  provide  useful  information,  if 
correctly  interpreted. 

3.  Scaling  the  geometric  factor  is  not  important  when  a 
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displacement  is  completely  dominated  by  viscous  forces 
(pseudo-stable  di splacement ) . 

4.  The  magnitude  of  the  linear  Leveret t  number  should  be  of 
the  order  of  0.01,  if  the  displacement  is  to  be 
stabilized.  Furthermore,  'stabilization'  of  the 
floodfront  may  be  defined  as  that  point  where  the 
recovery  is  within  5%  of  the  steady-state  recovery  as 
predicted  by  BucKley-Leveret t  theory. 

5.  The  interpretation  of  experimental  data,  obtained  for 
the  oil -wet  system  at  low  displacement  rates,  is 
complicated  by  the  effects  of  gravity  underride,  by 
frontal  instability,  and  by  a  non-stabi 1 ized 
displacement . 

6.  The  vertical  displacements,  conducted  at  low  superficial 
velocities,  revealed  that  the  behaviour  of  the 
waterflood  is  extremely  sensitive  to  any  heterogeneities 
in  the  porous  medium  and  to  the  'boundary  conditions'  of 
the  system. 

The  stability  of  the  displacements  was  also 

investigated.  The  following  conclusions  were  drawn. 

1.  Water-wet,  residual  water  saturation  system. 

a.  The  onset  of  instability,  as  predicted  by 
first-order  perturbations  theory,  at  Ns  =  13.56, 
(Peters,  1979)  was  verified.  Furthermore,  the  onset 
was  independent  of  the  geometric  factor. 

b.  Due  to  experimental  limitations,  it  was  not  possible 
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to  obtain  data  in  the  pseudo-stable  region,  where 
the  recovery  is  insensitive  to  viscous  fingering. 

c.  Displacement  tests  conducted  at  very  small  values  of 
the  stability  number  confirmed  that  the  recovery  was 
independent  of  this  group. 

2.  Oil -wet,  no  residual  water  saturation  system 

a.  The  recovery  at  the  onset  of  instability  was  seen  to 
be  different  for  the  various  sandpack  lengths 
considered.  This  was  attributed  to  the  effects  of 
gravity  underride. 

b.  The  onset  of  a  pseudo-stable  region  was  determined 
experimentally  to  occur  at  a  stability  number  of 
900. 

c.  The  recovery  is  independent  of  the  stability  number 
and  geometric  factor  in  the  pseudo-stable  region. 
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Recommendat ions 

The  following  recommendations  to  the  experimental 
design  of  the  equipment  are  suggested. 

1.  The  determination  of  water  arrival,  instead  of  water 
breakthrough,  should  be  attempted.  This  could  be  made 
possible  by  the  use  of  a  transparent  end-cap,  by  the  use 
of  radioactive  tracers,  or  by  the  use  of  microwaves. 

2.  A  hydraulic  ram  should  be  constructed  in  order  that  the 
sandpacks  may  be  extracted  from  the  long  coreholders. 
This  would  allow  the  capturing  of  viscous  fingers  on 
photographic  film. 

3.  The  displacing  and  displaced  fluids  should  be  of  'equal7 
densities  in  order  to  eliminate  gravity  underride. 

The  following  recommendations  for  further  study  are 
suggested . 

1.  Displacement  tests  should  be  conducted  in  rectangular 
sandpacks  in  order  to  verify  that  scaling  the  geometric 

factors  (  L  /L  ,  L  /L  )  is  not  necessary  for  stable  and 

x  y 

pseudo- stable  displacements. 

2.  A  set  of  scaling  groups  should  be  derived  by  the  method 
of  modified  inspectional  analysis  for  a  radial 
coordinate  system.  Verification  of  this  set  of  scaling 
groups,  with  a  reservoir,  should  be  undertaken. 

3.  The  stability  of  an  oil-wet  system,  containing  a 
residual  water  saturation,  should  be  investigated. 

4.  There  is  a  need  to  modify  the  stability  theory  in  order 
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to  be  able  to  predict  the  onset  of  the  pseudo-stable 
region.  It  is  suggested  that  the  investigation  parallel 
the  work  of  Bellman  and  Pennington  (1954). 

5.  A  series  of  tests  should  be  conducted  in  the 

pseudo-stable  displacement  region  for  a  water-wet 
system.  This  could  be  accomplished  either  by  increasing 
the  mobility  ratio,  or  by  decreasing  the  interfacial 
tension,  one  order  of  magnitude.  These  tests  should 
include  several  geometric  factors. 
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Appendix  A--Der ivat ions 

A.1  Derivation  of  Equation  2.16 

Restating  Equations  2.1  to  2.4  and  2.6 


w 


W 


w 


(grad  P  -  p  g) 
3  w  Mw 


v 


K  (5  ) 

- —  (grad  P  -  p  g) 

y^  o  o 


o 


as. 


div  v 

,  W 

■  *  at 

w 

as 

div  v 

0 

P  (s  )  = 
c  w 

p  -  p 

0  w 

Equations  A.  1.1  and  A. 1.2  may  be  rearranged  to: 


v 


w 


grad  Pw  -  Pw  9 


K  ($  )/y  =  grad  Po  "  po  g 

cr  w'  w 


Subtracting  Equation  A. 1.6  from  A. 1.7  results  in 


w 


TO  ■  'WTO  '  '  »”d  p-  *  -  v 


A.  1  .  1 

A.  1  .2 

A.  1  .3 

A.  1  .4 

A  .  1  .  5 

A.  1  .6 

A.  1  .7 

A.  1  .8 
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Ap  =  p  -  p 
M  MW  M0 


A.  1  .9 


and  reca 1 1 i ng  that 


grad  P  (Slf)  =  grad  -  grad  P 
c  w  3  o  3  w 


A.  1  .  10 


Equation  A. 1.8  becomes: 
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w17^ 


K  ( S,  ,)/u  =  grad  -  Ap9 


0V  W "  0 


A.  1  .  1  1 


Because  the  fluids  are  incompressible  and  immiscible  then: 


-y  -y  -> 

V  =  Vif  +  v 

W  0 


A.  1  .  12 


Using  this  relation,  Equation  A. 1.10  may  be  written  as 
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w'A'w  "oww"iJo 


K  (SJA 


=  grad  Pc(Sw}  -  Apg  +  K0(Sj/p‘0 


A.  1  .  13 


Defining 


s 


W'%  +  VswK' 


__  t,  +  MV 

\  Kw(Sw)  ^0  / 


-1 


A.  1  .  14 
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then 


t...  ■  us,  K-(S>’ 


W  W'  W'  y 


0 


KjS  ) 

qrad  P  (S  )  -  F  (S  )  — — —  Apg 
a  cv  w'  wv  w'  y 


A.  1  .  15 


Substituting  the  above  into  Equation  A. 1.3  results  in: 


K  (S  ) 

div  (  F  (S  )  — — —  grad  P(S,  1 

1  W  W  \in  c  w 


)  • d”  ( 


k  (s  ) 

F  (S  )  0  w  Apg 

Wv  Wy  yQ 


) 


+  div  (  Fw(Sw)v  )  +  *  3t 


w 


=  0 


A.  1  .  16 


From  vector  calculus 


div  lFw(sw)v)  =  FW(SW)  div  v  +  grad  Fw(Sw)  •  v 


A.  1  .  17 


Addition  of  Equations  A. 1.3  and  A. 1.4  and  noting  that: 


S  +  =  1 

w  o 


A.  1  .  18 


shows  that 


div  v  =  0 

Hence  Equation  A. 1.16  results  in: 
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A.  1  .  19 
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grad  PC(SW) 
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Applying  the  chain  rule  to  the  partial  derivatives  yields: 


grad  P„(SU1)  =  — ^  w  grad  Sw 


cv  w' 


w 


A. 1 .21 


and 


dF  (S  ) 

grad  F  (S  )  =  w  w 

3  wv  w' 


dS. 
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grad  S 
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A. 1 .22 


Then,  Equation  A. 1.20  may  be  written  as: 

K  (S  )  dP  (S  ) 


divl  Sw 


)  .  «,( 


Wv  W'  y 


Apg 


o 


dF  (s  )  as 

+  v  .  grad  Sw  +  *  ^  =  0 


dS 
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A.  1 .23 


Finally,  the  introduction  of  relative  and  base 
permeabilities  yields  the  desired  result: 


K  .  /  dP  (S  ) 

—  div  (  F  (S  )K  (S  )  — grad  S 
1  wv  w'  rov  w'  dS. .  ^  w 
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d1vlWW  Ap? 


-  dW  „  c  _  3Sw 

+  v  •  -J5- —  grad  Sw  +  *  ^ 
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A. 1 .24 
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A. 2  Derivations  of  Scaling  Groups 

The  following  section  presents  the  details  of  the 
derivation  of  the  scaling  groups  by  the  methods  mentioned  in 
Chapter  2. 

A. 2.1  Dimensional  Analysis 

Loomis  and  Crowell  (1964)  and  Greenkorn  (1964) 
presented  a  dimensional -analysis  approach  to  determine  the 
scaling  groups  for  immiscible  flow,  and  miscible  flow  with 
heat  transfer,  respectively.  The  former,  however,  omitted 
time  in  the  analysis. 

Presented  below  is  a  derivation  of  the  scaling  groups 
for  the  immiscible  fluid  displacement  problem.  The 
derivation  follows  closely  the  report  of  Loomis  and  Crowell 
(1964),  except  for  a  few  small  modifications  tailored  to  the 
definitions  of  capillary  pressure  and  relative  permeability 
adopted  in  this  study. 

The  variables  which  control  the  displacement  process 

are : 

1 )  Geometric: 

a)  length  of  the  system 

b)  thickness  of  the  system  or  diameter 

c)  angle  of  dip 

2)  Sand  properties: 

a)  permeability  to  oil  at  residual  water 
saturation 


' 

. 


- 
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b)  permeability  to  water  at  residual  oil 
saturation 

c)  porosity 

3)  Fluid  Properties: 

a)  water  viscosity 

b )  oi 1  vi scosi ty 

c)  water  density 

d)  oil  density 

e)  capillary  pressure  normalizing  parameter 

4)  Others: 

a )  time 

b)  gravity 

c)  superficial  velocity 

Upon  neglecting  inertial  forces,  the  variables 
p  ,  p  ,  and  g  occur  only  in  the  combination  Apg ,  where: 

AP  =  pw  -  p0  A. 2.1 

The  functional  relationship  of  these  variables  may  be 
written  in  the  following  implicit  form: 

f  (L,  h,  a,  KQr,  K^,  yw»  yQ»  Apg,  cr,  V,  t)  -  0  A. 2. 2 

Applying  the  principle  of  dimensional  homogeneity  results 

in: 

mo*oto  i  Lai  ha2  Oo  k^3  ,0  ua5  y»6  (ipg)a7  0as  va9  ta10  A.2.3 
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where  a: ,  a2 ,  a3,  .  a  1Q  are  arbi trary  constants . 
Introducing  the  dimensions  of  the  variables  into  the  above 
equation  yields: 

m0£0t°  4  £a!  Za 2  42a3  nZab  (m  t"1  )**  (m  d"1  t"1)36 

(m  d'2  t"2)  (m  d'1  t"2)  8  (d  t"1)  9  tai°  A. 2. 4 

-  u|(^5"*’^6'*’  a7  +  ag  )  ^  (al  +  a2  +  2a  3  +  26^-  95”  8  5  “237-33+89) 


t 


("  a5  "  a6  “  2a7  -  3g  +  310) 


A. 2. 5 


Equating  co-efficients  of  m,  &  ,  and  t  results  in: 


a5  +  a6  +  a7  +  a8  “  0  A  .  2 . 6 

3l  +  82  +  283  +  234  -  85  -  83  -  2a7  -  83  +  8g  =  0  A. 2. 7 

-  85  -  8g  -  2a7  -  8g  +  Sjo  =  0  A. 2. 8 

It  must  be  established  that  the  above  set  of  equations  is 
independent  in  order  to  determine  the  minimum  number  of 
pi-terms.  In  matrix  form,  Equations  A. 2. 6  to  A. 2. 8  may  be 
written  as: 


* 


' 


0  0  0  0 
112  2 
0  0  0  0 


1  10  0 
-2-110 
-2  -2  -1  1 


A. 2. 9 
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1  1 

-1  -1 

-1  -1 


p  dl 

a2 

• 

it 

r  o 

0 

• 

• 

.  aioJ 

i'ol 

A  theorem  from  matrix  algebra  ( Langhaar ,  1951)  states  that 
the  number  of  independent  equations  will  be  equal  to  the 
order  of  the  highest  order  determinant  of  the  matrix  which 
is  different  from  zero.  The  value  of  the  determinant  may  be 
positive  or  negative. 

The  highest  order  of  the  determinant  for  the  above 
co-efficient  matrix  is  3.  It  may  be  shown  that  it  is 
possible  to  write  the  3-d  determinant  which  has  a  non-zero 
value.  Thus,  from  Buckingham's  (1915)  pi -theorem,  the 
minimum  number  of  pi -terms  is  (the  number  of  variables,  12, 
less  the  number  of  independent  equations,  3)  9. 

A  possible  solution  to  the  above  system  of  equations, 
expressing  three  of  the  co-efficients  in  terms  of  the 
others,  is: 


a i  =  2a2  ~  ^a3  “  2ai+  +  a7  -  ag 


A. 2. 10 


a 5  =  -  a6  -  a7  -  a8 


A. 2. 1 1 


a1Q  =  a7  +  a8  +  ag 


A. 2. 12 


Substituting  into  Equation  A. 2. 3  results  in: 


- 
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m° £° t 0  £  i_(-a2-233-2a4+a7-a9)  .a2  o  i<a3  i/a4  4o 

or  wr  ^ 

u^-a6-a7-a8)  ya6  (Apg)a7  aa8 

ya 9  ^.(a7+a8+a9) 


which  upon  rearranging  yields: 

d  (L_1h)a2  (L'2K  )a3  (L'2K  )a-  )a6 

v  '  v  or  wr7  K  vj  o 


(Ly^Apgt)9?  (L_1Vt)a9  a°  *<> 


From  Equation  A. 2. 14  arise  the  following  pi-terms: 


h 

"l  =  t 


tt2 


^3 


\jr 

L2 


7T4 


"5 


LApgt 
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"6 


.2.13 

.2.14 

.2.15 

.2.16 

.2.17 

.2.18 

.2.19 


A. 2. 20 
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A. 2. 21 


tt8  -  a 


A. 2. 22 


tt9  =  <j>  A.  2. 23 

The  foregoing  pi -terms  may  be  brought  to  the  form  below 
because  it  is  permissible  to  multiply,  divide  or  invert 
these  terms,  provided  that  a  complete  set  of  pi -terms 
contains  every  variable  at  least  once,  thus: 


"l 
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or 


L2 


*3 
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wr 


L2 


7T4 


Kwr  yo 
^or  uw 
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K  Apg 
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Kwr« 


IT? 
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tt8  =  « 


TTg  -  4> 


A. 2. 24 

A. 2. 25 

A. 2. 26 

A. 2. 27 

A. 2. 28 

A. 2. 29 

A. 2. 30 

A. 2. 31 

A. 2. 32 
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A. 2. 2  Inspectorial  Analysis 

The  equations  describing  the  fluid  displacement 
problem,  as  given  by  Equations  2.1  to  2.6,  may  be  combined 
(see  Appendix  A.1)  to  yield: 

K  .  r  dP  (S  )  1 

-  div  F  (S  )  K  (S  )  — 5?-  -w--  grad  S 

u_  wv  w7  rov  w7  aS  y  w 

o  |_  w 

-  —  div  I" F  (s  )  K  (S  )  Aog)l 

yQ  ^  wv  w7  r o  w  y  J 

dF  (S  )  9S 

+  7  •  -V1-  9radSw+  *3T  =  5  A'2’33 

w 

Upon  substituting  the  following  capillary  pressure 
model  (Bentsen  and  Anli,  1976); 


P  (S  )  = 

cv  w7 


air  (S*) 
c 


+  p 


A. 2. 34 


and  using  the  end-point  relative  permeability  definitions  as 
a  function  of  normalized  water  saturation: 


K  (S*)  = 

rvr  7 

K  (S*)  = 

rcr  7 


Vs*> 

V 

K  (s*) 

0V 

K 


A. 2. 35 


A. 2. 36 
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Equation  A. 2. 3 3  becomes: 


crK 


or 


di  v 


dTT  (S*) 

Fw(S*>  Kro(S*>  -asr-  9rad  s* 


K 


or 


o 


di  v 


Fw(S*}  Kro(S*}  Ap5 


dFw(s*) 


8S* 


+  V  ‘  ~d$*  9rad  s*  +  <J>  TF  =  0 


A. 2. 37 


Expanding  Equation  A. 2. 37  in  cylindrical  coordinates  (r ,x, b) 
and  retaining  the  gravity  term  results  in: 


ctK 


or  \  ]_  _8_ 
I  r  8r 


d*r(S*)  as* 

r  F  (S*)  K  (S*)  — -  —  - 
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wv  ‘  rov  '  dS*  83 


+  ± 
8X 


<MS*)  as* 

F  (S*)  K  fS*)  _ - _ 

V*  '  ro^  '  dS*  8 x 


ApgK 


u. 


or  J  I  J_ 
r  8r 


r  Fw(S*}  Kro(S*> 


]_  _3_ 

vT. 


F  (S*)  K  (S*) 
w  ro '  ' 


3X 


Fw($*)  Kro(S*) 


+ 

*  3t 


v  IS!  +  \  +  v  |§i 

r  8r  r  83  x  8x 


dFw(S*) 

dS* 


=  0 


A. 2. 38 


The  next  step  in  an  inspectional  analysis  is  to  define 
certain  dimensionless  variables.  After  the  manner  of  Perkins 
and  Collins  (1 960 ) : 


107 


T 


|  V|  t 

*L(l  -  S 


-  S  .) 

or  wi ' 


and 


A. 2. 39 

A. 2. 40 

A. 2. 41 

A. 2. 42 


S*  = 


S  -  S  . 
w  wi 


1  -  S  -  S  . 

or  wi 


A. 2. 43 


Substitution  of  Equations  A. 2. 39  to  A. 2. 43  into  Equation 
A. 2. 38  yields: 


a  K 


or  )  1 


o  (  a2\p 


+ 


ApK 


8 

dip 

- 

1 

a2\p 

2 

1 

3 

L2 

3C 

1 

8 

aip 

8^ 

1 

8 

L 

8^ 

diT  (S*) 

Kro^  35*“" 


dip 


8 

36 


Fv/S*>  Kro^S*)  -Is* 


d\  (S*}  3S* 


36 


<MS*) 

F  (S*)  K  (S*)  — 
wv  '  ro  ‘  dS* 


3? 


*F,j(S*)  Kro(S*)  | g |  cosa 


Fw(S*)  Kro<S*)  sina 


3$*  +  lV^f6  1  3S*  , 


Ivlf 


dip 


a  ^36 


I  3S* 

3? 


dFw(S*) 

dS* 


|vl  8S* 
L  9x 


=  0 


A. 2. 44 
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The  above  equation  may  be  rearranged  as: 


+ 


5  3? 


dF  (S*) 
wv  ‘ 

dS* 


+ 


9  S* 
9t 


0 


A. 2. 45 


From  the  dimensionless  equation  above,  the 
dimensionless  groups  which  arise  are: 

1)  Independent  dimensionless  variables 


B »  t 


2)  Dependent  dimensionless  variables 


S*, 


. 
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3)  Independent  constant  or  similarity  groups 


ctK 


wr 


a  K 


wr 


K  u  K  Apg 
or  w  wr 


|v|uwD 


I  . 

v  y  L 
1  Mw 


D  '  'Vo 


V  U 


w 


KJS*>>  Kro(S*>’  \<S*> 


A  similar  analysis,  using  the  Leverett  J-function  to 
represent  the  capillary  pressure,  results  in  the  following 
simi 1 ar i ty  groups : 


L 

D 


*^or  yw 
*^wr  uo 


Kwrip9 


J 


KJS*),  Kro(S*),  J(S*) 


/ 


no 


A. 2. 3  Modified  Inspect ional  Analysis 

Restating  Equations  2.20  and  2.21 


t 


W  Kob 


w 


I  ~ ^ 

V 


K  (S.,)  C  W 


rcr  w' 


dPc(Sw) 

3s —  grad  s 


W 


w 


F  (S  ) 
wv  w' 


K 


ob 


v 


1  w  1 

rov  w7 

W  v  W '  |  •>  | 

1  vl 

0 

1  vl 

di  v  f 


w 


=  _  JL 

Iv1 


A. 2. 46 


A. 2. 47 


To  complete  the  description  of  the  fluid  displacement 
problem,  an  appropriate  set  of  initial  and  boundary 
conditions  is  required. 


w 


w 


=  S 


Wl 

=  1  -  S 


or 


-> 

f 

f 


w' 


w1 


=  0 
=  1 


1) 

2) 

3)  0  <  Sw  <  1  f?J  =  1,  t  =  0 

Condition  1  states  that  initially,  throughout  the  sandpacK, 
a  residual  water  saturation  equal  to  SW1-  exists;  condition 
2  implies  that  the  porous  medium  is  flooded  out  to  a 
residual  oil  saturation;  and  finally,  condition  3  represents 
a  '  jump  discontinuity'  in  the  fractional  flow  at  the  inlet. 

This  completes  the  description  of  the  problem.  For  each 
variable  in  the  problem,  a  dimensionless  variable  may  be 
constructed  by  dividing  the  original  variable  by  some 
arbitrary  reference  quantity  of  the  same  dimensions,  viz: 


c  E  3T 

xo 


A. 2. 48 


e 


T  = 


S* 
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A. 2. 49 

A. 2. 50 

A. 2. 51 

A. 2. 52 


where  x  ,  r  ,  3  ,  tQ,  ,  and  Bq  are  arbitrary  reference 
quantities.  Equations  A . 2 . 48  to  A. 2. 52  are  then  substituted 
into  the  fluid  displacement  equations,  initial  and  boundary 
conditions.  Consider  first  the  conditions: 


S*  Bo  +  Ao 


A. 2. 53 


or 


S* 


f  I  =  0  A. 2. 54 

w1 


Also, 


S 

w 


S*  B  +  k 
0  0 


A. 2. 55 


Hence 


-V 


S* 


A. 2. 56 
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Now  for  | f w |  -  0,  S*  -  0  ;  and  from  Equation  A. 2. 54 

A  =  S  .  A. 2. 57 

o  wi 

Also  for  | f  |  =1,  S*  =  1  ;  and  from  Equation  A. 2. 56 

w 


B 


o 


A.  2. 58 


The  same  may  be  done  for  Equations  A. 2. 46  and  A. 2. 47; 
but,  before  proceeding,  a  dimensionless  form  of  the  gradient 
and  divergence  must  be  examined.  Recall  that  in  a 
cylindrical  coordinate  system: 


grad  h  e 


36 


A. 2. 59 


and 


div  t 


1 

r 


+ 


30  £0 


+ 


r 


_3_ 

3X 


A. 2. 60 


Upon  substituting  for  the  variables  defined  by  Equations 
A. 2. 48  to  A. 2. 50,  Equations  A. 2. 59  and  A. 2. 60  become: 


and 


grad  h  =  — 


_J _ 3_  -  +  J 

w  3e  e  xo 


d  ,  I  L. 

- —  E  +  -  -  h 


A. 2. 61 


di  v  £  = 


r  $ 
o 


r  $ 

+  -2-  a 
x0  3C  C 


A. 2. 62 


. 
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Recall  Equation  2.24: 


pc(S*)  =  <rirc(S*)  +  Pd 


A. 2. 63 


to  be  used  below. 

Equations  A. 2. 46  and  A. 2  47  may  now  be  written  in  the 
following  form: 


-y 

f 


W 


Fw(S*>  Kor  <MS*) 

^  V7- Kro^  ° -hr~ 


1  aS*  - 

r  a$  $ 

o 


1  as*  -  L  1  as*  - 

+  — _ -  £  +  tr-  T7=—  C 


r  B  $  3e 

0  0 


x0  3? 


F  (S*)  K 

w  Ap9 +  Fv/S*> 


-y 

V 


I 

V 


I 

V 


A. 2. 64 


and 


r  $ 
o 


'  i  a f  r  $  af  ' 

l(*f  )  +  J - +  -2 _ «£ 

a$  v  w$y  6q  ae  XQ  a? 


<J)B 


o  as* 


i^i t0  3T 


A. 2. 65 


Rearranging  the  above  equations  leads  to: 


+  dTT  (S*) 

f  =  F  (S*)  K  (S*)  --4;-*— 

w  w  ro  dS* 


aK 


or  aS*  * 

$ 


i->-i  a$ 

v  hi  r 
1  o  o 


aK 


or  1  as*  - 

e  + 


aK 


i  -*■  i  n  $  a  e 
v  u  r  $ 

1  1  Mo  o  o 

-y  K  Apq 

1 - °i— Ik  (s*) 

-y,  |  ■>  I  ro v 

V  V  hi_ 


or  aS*  £ 


->i  a  r 

V  hi  X  ^ 
1  0  0 


F  (S*) 


A. 2. 66 
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and 


I  v|  t 


4>B  r 
oo 


—  —  ($f  ) 

$  3$  v  W$' 


V  t 


0 


cj)B  r  g 
*  o  oo 


1 

$ 


3f. 


We 


3e 


vlt 


4>B^x 
r  0  0 


3  f. 


wc 


3C 


as * 
aT 


A. 2. 67 


In  the  above  two  equations,  the  parameters  and 
reference  quantities  have  been  grouped  together  to  yield  a 
description  of  the  problem  in  terms  of  dimensionless 
variables.  These  dimensionless  groups  may  then  be  equated  to 
unity  to  form  a  set  of  algebraic  equations.  These  algebraic 
equations  may  then  be  solved  to  determine  the  reference 
quantities,  viz: 


and 


aK 


or 


=  1 


vIVo 


aK 


or 


v  u  0  r 
'  Mo  o  o 


=  1 


aK 


or 


1*1  Vo 


=  1 


hence  rQ 


hence  80 


hence  xq 


=  1 


A. 2. 68 

A. 2. 69 

A. 2. 70 


v  t 


o  


4>B  r 
oo 


=  1 


hence  t 


d>(l-S  -  S  .)  aK 
or  wi '  or 


o 


A. 2. 71 


v  y 


o 


Although  identical  expressions  were  found 

for  x  and  rn,  this  is  not  the  case  in  general.  These 

o  u 

quantities  should  not  be  forced  to  be  identical  as  this 
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would  lead  to  a  loss  of  valuable  information. 

Now,  for  the  sake  of  convenience,  the  following 
dimensionless  groups  are  defined: 


and 


A. 2. 72 


A. 2. 73 


where  a  is  the  core  radius  and  L  is  the  length  of  the 
sandpack.  Recall  Equation  2.11  in  terms  of  the  end-point 
permeabi 1 i t ies : 


K 


or 


y 


w 


A. 2. 74 


The  dimensionless  variables  defined  by  Equations  A. 2. 48  to 
A . 2 . 52  now  become : 


$ 


r_ 

a 


e 


x 

L 


Ivlt 


♦l1  "  Sor  “  Swi^ 


A. 2. 75 

A. 2. 76 


T 


A. 2. 77 
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or 


T 


*0'  -  s 


A. 2. 78 


The  physical 
is  not  an  obvious 


interpretation  of  the  variables  $  ,c  and  T 
one.  Defining  the  following  variables: 


and 


A. 2. 79 


_  x  .  cNc£  A.  2. 80 

L  M 

r 


t  -¥  | 

V 


TN 


C£ 


<J>L(1  -  S 


or 


-  S  .) 
wi ' 


A. 2. 81 


It  is  now  evident  that  V  ,  K  and  x  have  the  physical 
interpretation  of  dimensionless  radius,  dimensionless 
length,  and  dimensionless  time,  respectively.  Furthermore, 
dimensionless  time  corresponds  to  pore  volumes  of  fluid 
injected  into  the  system.  As  a  consequence  of  the 
introduction  of  Ncr  and  Nc£  ,  the  number  of  parameters  is  no 
longer  a  minimum.  That  is,  N  and  Nc£  are  redundant. 

The  above  derivation  of  the  reference  quantities  is 
based  on  a  Lagrangian  formulation  of  the  initial  and 
boundary  conditions.  The  result  is  to  yield  a  differential 
equation  which  does  not  account  for  the  length  and  diameter 
of  the  system.  These  dimensions  are  accounted  for  in  the 


' 
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non - norma  1 i zed  dimensionless  lengths, 


and 


0  <  c  < 


1 


0  <  $  < 


1 


and  $  : 


The  derivation  may  also  be  carried  out  using  the 
Eulerian  form  of  the  initial  and  boundary  conditions.  Two 
possible  definitions  will  arise  for  the  arbitrary  reference 
quantities,  viz: 


1) 


The  first  definitions  are  the  same  as  the  ones  which  arise 
from  the  Lagrangian  formulation.  As  a  result  of  using  the 
second  definitions,  the  differential  equation  will  account 
for  the  dimensions  of  the  porous  medium  and  the 
dimensionless  distances  will  be  normalized,  viz: 


0  <  if,  <  1 


and 


0  <  5  <  1 


As  a  consequence  of  the  introduction  of  the  definitions  for 


• 
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^  and  K  (Equations  A. 2. 79  and  A. 2. 80),  essentially 
equivalent  results  may  be  obtained. 

The  fractional  flow  equation  and  the  equation  of 
continuity  may  now  be  written  in  terms  of  ip  ,  5  and 


and 


,  viz: 

-> 

f  = 
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A. 2. 82 
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C&  3t 


A. 2. 83 


Equations  A. 2. 82  and  A. 2. 83  may  be  simplified  to: 


fw  =  G(S*)  -  NCJ,C(S*) 


2L  / 3S*  «  ,  1  3S*  s  \  j.  3S* 
'D-(3F*  +  ?3S“S)+  3? 


A. 2. 84 
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A. 2. 85 


where 
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N  1 
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r  j 

r*pg 
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A. 2. 86 


A. 2. 87 
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dir  (S*) 

C(s*)  =  -hr- 


s*  = 


s  -  s  . 

w  wi 


-  S  -  S  . 
or  wi 
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*^wr  yo 
*Sdt  uw 


Fv/S*> 
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M  K  (S*)  +  K  (S*)' 
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A. 2. 88 

A. 2. 89 


A. 2. 90 


A. 2. 91 
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Appendix  B- -Exper imenta 1  Results  From  the  Literature 


TABLE:  B.1 

SUMMARY  OF  SANDPACK  PROPERTIES 
DATA  FROM  PETERS  (1979) 

WATER-WET  SYSTEM 


RUN 

LENGTH 

DIAMETER 

POROSITY 

ABS.  PERM. 

PERM.  @>SWI 

SWI 

(CM) 

(CM) 

(%) 

(DARCY) 

(DARCY) 

(%) 

1 

22.9 

4.809 

38.98 

16.43 

15.82 

1  1  . 

2 

23.6 

4.844 

36.37 

18.36 

17.99 

8. 

3 

112.8 

4.961 

38.38 

21.89 

18.05 

9. 

4 

23.7 

4.809 

38.94 

18.93 

18.33 

1  1  . 

5 

110.5 

4.825 

38.91 

20.52 

18.50 

10. 

6 

110.0 

4.805 

37.61 

21  .91 

19.40 

10. 

7 

23.6 

4.844 

37.49 

14.  19 

11.03 

15. 

8 

116.1 

4.961 

34.66 

18.28 

15.57 

8. 

9 

113.0 

4.958 

39.62 

22.99 

19.90 

8. 

10 

115.9 

4.958 

35.86 

19.22 

15.49 

12. 

1  1 

22.8 

4.843 

37.52 

16.23 

15.  19 

1  1  . 

12 

110.4 

4.805 

37.46 

20.85 

18.48 

10. 

13 

115.9 

4.972 

37.80 

18.54 

15.53 

9. 

14 

110.0 

4.805 

38.11 

22.50 

18.62 

8. 

15 

1  12.8 

4.972 

35.48 

22.72 

20.90 

9. 

OIL 

-WET  SYSTEM 

RUN 

LENGTH 

DIAMETER 

POROSITY 

ABS.  PERM 

(CM) 

(CM) 

(%) 

(DARCY) 

19 

22.9 

4.809 

41.33 

17.26 

20 

23.6 

4.844 

37.89 

15.05 

21 

22.9 

4.809 

39.64 

18.02 

22 

110.5 

4.825 

38.46 

22.63 

23 

23.5 

4.843 

38.25 

16.20 

24 

23.5 

4.809 

38.04 

16.34 

25 

115.6 

4.972 

37.01 

18.36 

26 

23.7 

4.809 

38.20 

15.62 

27 

22.9 

4.844 

39.82 

16.  18 

28 

110.0 

4.825 

38.  13 

24.79 

29 

115.7 

4.961 

37.49 

18.62 

2 

7 

2 

3 

8 

5 

3 

3 

6 

5 

2 

0 

3 

9 

9 


• 

TABLE:  B.2 

SUMMARY  OF  DISPLACEMENT  TESTS 
DATA  FROM  PETERS  (1979) 
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Appendix  C--Sample  Calculation  of  Extrapolated  Data 

This  section  demonstrates  the  calculations  for  Point  E  of 
Figure  4.1  for  the  inverse  linear  Leverett  correlation.  The 
magnitude  of  the  stability  number,  Ns,  corresponding  to  the 
point  where  the  viscous  forces  completely  dominate  the 
displacement  process,  is  900.  Thus, 


Solving  for  the  velocity  yields: 


V  = 


900C*yK 


(H 


C .  2 


u  D2 


Recalling  the  definition  of  the  inverse  linear  Leverett 
number : 


1 


V%L 


Lt 


i  y  (  K<J> ) 


0.5 


C .  3 


and,  upon  substituting  Equation  C.2,  gives  the  result: 


Lt, 


900C*YK 


y  L 
w 


Jw  / 

Equation  C.4  may  be  simplified  to: 


u, 


Y 


(K*) 


0.5 


C.4 


1  .  900C*  ( K/  ij> ) 0  ~ 


Lt, 


(H 


D2 


C .  5 


The  inverse  linear  Leverett  number,  in  a  water-wet, 
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connate-water-bearing  system  may  now  be  estimated  from 
Equation  C.5.  In  particular,  for  Point  E: 

C*  =  306.25 


K  =18  Darcys  =  1.78x10-7  cm2 

4>  =  0.36 

L  =  22.6  cm 

D  =  4 . 8 1 5  cm 

where  L  and  D  are  the  dimensions  of  a  sandpack,  and  <J>  and  K 
are  the  estimated  average  sandpack  properties.  With  these 
values  for  the  variables,  the  inverse  linear  Leverett  number 
is  1.86  x  10°.  The  same  type  of  calculation  may  be  carried  out 
for  any  other  case  desired. 


